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ABSTRACT 


In this study, shape preserving data driven rational cubic schemes are 
developed. A rational cubic piecewise function (quadratic denominator and 
cubic numerator) with two parameters was transformed to C! rational cubic 
piecewise function. Constraints were derived on free parameters by means 
of some mathematical derivations to train and demonstrate convex curve. 
The scheme, then, was advanced to partially blended rational bi-cubic 
function with eight free parameters which were controlled to ascertain 
convex surface. A numerical comparison with certain existing schemes 
manifested that the proposed method was economical. The proposed 
scheme was put into visualization of convex 2D and 3D data using 
MATLAB software packet. Additionally, the suggested approach produced 
a more visually appealing interpolating curve for scientific visualization for 
specific data sets. 


Keywords: convex data, convex surface, rational functions, shape 
preservation 


1. INTRODUCTION 


A considerable amount of raw data that represents the experiment's 
outcomes is typically obtained by the scientists when the experiments are 
carried out in the lab. For instance, the data in a chemical experiment might 
be the distillation curve, whereas in a biology experiment, it could be the 
movement of the relevant cell within the system. Additionally, any 
experiment may have error or noise in the data collecting procedures. The 
type of error could be systematic, biased, or random. Other methods of data 
processing, such as data fitting and data smoothing or filtering were applied 
to these types of data. Interpolation is a preferable alternative for the data if 
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the data is clean, which means there are probably no errors. However, if the 
data contains noise, approximation, smoothing, or de-noising may be 
advantageous to gain fruitful results. There are several popular techniques 
for processing the data including linear regression, polynomial regression 
(or fitting), any non-polynomial fitting, cubic splines, and smoothing 
splines. Interpolation, on the other hand, addresses the problem of creating 
curve or surface that must also maintain the data's geometric shape when it 
is approximated which is subject of the current research. 


Science and engineering, both rely on these continuous visual 
representations of discrete scientific data in particular which are coherent 
with intrinsic properties of the data. Computational methods used in these 
realms are terms as shape preserving algorithms. The design of shape 
preserving data embedded curves and surfaces lies at the heart of several 
scientific disciplines in the digital age including computer-aided geometric 
design (CAGD), computer graphics, and computer-aided design (CAD). 
These applied areas combine via computational software packets, ideas 
from numerical techniques, visualization, differential geometry, and linear 
algebra. Surface and curve drawing schemes offered architects a 
competitive edge over traditional drafting tool and optimized the likelihood 
of realizing their ideas by exhibiting their thoughts on computer displays in 
remarkable short amount of time. 


Moreover, the design of machine software, robotics, and chemical 
processes are among the many fields in which shape-preserving approaches 
are widely used. For designing ceramics, cutting wood for wooden furniture 
and creating ornaments, a lathe is utilized. The design of the lathe path is 
made by using shape-preserving software. Shape-preserving interpolation 
techniques are employed in both, chemistry and robotics, to plan the path of 
the robots. In chemistry, they are used for PARFAC modelling of the 
fluorescence data. There are more applications in the areas of signal 
processing and photogrammetric data manipulation. Convexity, 
monotonicity, and positivity are the most dominant data types. All of these 
applications have to be retained in the process of visualization. 


The usual data interpolation algorithms present the data smoothly, 
however, fail to intact its shapes. The existing shape-preserving algorithms 
restrain the derivatives at the knots to envisage convex, monotone, and 
positive data as convex, monotone, and positive surfaces and curves and 
operate effectively with specific conditions, for instance, the solution of 
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ordinary and partial differential equations cannot be interpolated by these 
shape-preserving methods. Another strategy is to add additional knots 
between the existing data points. Even though these techniques maintain the 
forms of the data, the amount of computing work has grown larger. 


The motivation behind the current research was to devise an efficient 
and economic algorithm that not only delivers smooth visualizations but 
confirms convexity in the visual attribute. This research addressed the 
convex surface and curve interpolation which has become a fundamental 
element of numerous applications in several areas, such as airplane and 
automotive industries, architecture, mechanical engineering, and computer 
graphics. Likewise, convexity in yield curve measures the extent at which 
price of a bond changes with the rate of interest prevalent in market. It 
benefits the stakeholders to anticipate price of a bond if market interest rate 
changes. Furthermore, indifference curve in economics which ranks a 
consumer’s preference for a combination of goods represents convex 
structure. 


Shape conserving convex interpolation was well explored in past few 
decades and some eminent contributions are [1—26]. The control points of 
the boundary of 2D shape are optimized using B’ezier curve in [1] and [2]. 
BFGS diagonal approximation method was employed in [3] to solve B 
spline curve fitting problem. Surface reconstruction, using B spline and 
BFGS diagonal approximation, was discussed in [4]. Kvasov [6] presented 
two algorithms to preserve the monotonicity and convexity of the discrete 
data by using weighted C! quadratic splines. Verlan [7] proposed an 
algorithm of C? interpolation utilizing splines of arbitrary degree for 
discrete data to achieve the convexity conservation of data. Moreover, in 
[8—10] convexity preserving techniques were global. 


Brodlie and Butt [11] presented a piecewise cubic C! interpolant to 
preserve the convexity of the data. Hussain et al. [12] suggested a local 
shape preserving interpolation scheme through СС? bi-quadratic 
trigonometric spline. The scheme comprised of four free parameters in each 
rectangular domain. These free parameters were controlled so that no 
possible fluctuations could be seen. This scheme maintained the smoothness 
of positive and monotone 3D data. Sarfraz et al. [14] constructed a rational 
cubic function with two parameters in the vision of captivating positive 
data. Simple constraints were manipulated on two parameters in the 
illustration of rational cubic function and extended rational bi-cubic 
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partially blended function. The scheme was local with an approximation 
order of O(h?). Sarfraz et al. [15] developed the smooth shape preserving 
methods deploying piecewise rational cubic and the rational bi-cubic 
function to conserve the convexity for 2D and 3D data. The data dependent 
constraints on free parameters were evaluated to achieve the shape of 
convex data. The convexity conserving schemes yielded adequate degrees 
of freedom to facilitate the researcher and user. The convex shape of the 
data was not preserved by the techniques in [17] when the data was 
presented with derivatives. 


To retain the shape of convex data, Costantini and Fontanella [22] 
created a semi-global approach. However, this scheme was not economical 
as degree of interpolant in some rectangular patches became excessively 
high and the polynomial patches tended to be linear in x and/or y. The 
resulting surfaces, also, failed to attain a visually pleasant look. In order to 
create convex surfaces, Asaturyan [23] each rectangular grid was divided 
into nine smaller rectangles. This scheme was not local which means that 
by altering the data for one sub rectangle's edge in the x direction, a 
modification was made to the grid for all of the sub rectangle's edges that 
were originally situated in the x direction. The convex data visualization 
techniques created in [24] were global. 


In [25], the proposed approach kept the convexity of the surface along 
the grid lines, but failed to preserve the convexity inside the grids, leading 
to unfavorable flat areas because second order mixed partial derivatives 
vanished. A univariate C? cubic spline was utilized to preserve convexity. 
A C? quartic \ quadratic interpolant was used by authors in [26] to visualize 
monotone, convex, and positive data set. 


The convex shape of the data was not preserved by the techniques in 
[17, 25] when the data was presented with derivatives. 


The main focus of the current paper was to derive data dependent 
constraints on parameters in the description of piece wise rational cubic and 
bi-cubic functions in [14] to hold convexity. The research conducted for this 
purpose was discussed in several sections. Section 2 and Section 3 presents 
a review of piece wise rational cubic and rational bi-cubic function [14] to 
be employed. Convexity preserving constraints for curve and surface data 
were developed in Section 4. Section 5 provides implementation of rational 
scheme developed on 2D and 3D data sets and corresponding numerical 
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results were demonstrated in Section 6. Proposed research algorithm was 
concluded in Section 7. 


2. RATIONAL CUBIC FUNCTION 


The rational cubic function introduced by Sarfraz et al. [14] is explained 
in this section. 

Suppose that (uj, v), where j = 0,1,2, ..., n are the set of given data 
points defined over the interval [a, b] such that a = ug < uy, <и, <= < 
Un = b. The piecewise rational cubic function with two parameters 0; and 
o; defined over the each subinterval J; —[uj,,uj where j= 
0,1,2,3, ...,n — 1 is given by: 
Ти) 

LB 
0“ 
where 
PG) = ojvj (1 — 2? + lojv; + hjojdj + 2v,}(1 — £)2£ + {оу — 
һо)”  2vj41)(1 — ë) С “ауы, 
QjG) = 0301 — £ +201 — 2€ + of, 


_ (u-uj) 
= 


(2.1) 


y hj = Uj+1 — uj. 


The rational cubic function defined in Eq. (2. 1) is C! if the following 
characteristics hold, 


T,(uj) = v Тушаал) = Yar T (uj) = do Tejas) = du. 0.2) 
where ТО) signifies derivative with respect to ‘u’ and d; symbolizes the 
derivative values at the knots иу. Interestingly, taking o; = o; = 1, rational 
cubic spline takes the form of standard cubic Hermite spline in each 
interval J; = [uj, uj+1]. 

3. RATIONAL BI-CUBIC FUNCTION 


In this section, the rational cubic function in Eq. (2. 1) is expanded to a 
rational bi-cubic function specified over three-dimensional set of data 
points {(uj, vj, G): j = 0,1,2,...,n; | = 0,1,2, ..., m ]with the rectangular 
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grid М = [ug, Un] X |”, Vm]. Let i:a = ug < u, < u; €: < Un = b be 
the partition of [a,b] and #:c = Vo < vi < V3 < = < Vp = d be the 
partition of [c,d]. The rational bi-cubic partially blended function can be 
defined over each rectangular patch Jj = [u;j, uj+1] х [Vi vi44] as: 


T(u,v) = —EGNT, (3.1) 
where 

0 T (u, vj) T (u, vi+1) 
G=| Т(ш,т) T (uj, vi) T(v,visi) |, (3.2) 


T(u.,v) Tjon) Т(шаы, Vira) 

E=[-1 т0(6) nG)] 
N —[-1 2m) 41671. 
where 
To = (1- 7? (1 + 22), т, = 22(3 – 22), 4920-0 + 20), А = 
n^ – 2n). 

Furthermore, 4 = 9 ÜSISI, n= ee 0x 
n < 1, 
hy = щу — Uj, hi = Viya — 0. 

The rational cubic function T (u, vj), T (u, Vi41),T (uj, v), T (uy, a, v) 


similar to Eq. (2.1) are attained and specified over the boundary of 
rectangular portion as: 


Daag 0M 


Tw) - aD C 


(3.3) 


where 

Mo = 016, 

М, = Qj 1G), + hjg Gi, + 26, 

M, = ЖЕГЕН = hyo б лн + 2Gj41 Мз = ЕГЕТ 
0.0) = 03101 2)? + (1 2) + оё? 


with 0/1 and оу shape parameters. 
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Likewise, 

3 1 — 3—1 ІМ 
талда) = ERU T (34) 
where 


No = Она біз» № = 0j 141 білі + оза Са + 26; 1415 
№ = біз брі - Һауа GA ua + 26; із Ns = біле бу» 
Q2(€) = Qj 14101 —  )2 + (1-2) + Oj 14187 
with shape parameter Qj 141 and o; 144 . 
Further, we have, 


_ Xj 7 070]; 


T (uj, v) = р ° (3.5) 
where 
Jo = ёб h = булбул + h; Oj iG + 26, 
Ja = баби — AG Ghia. + 2G м. з = 851654; 
Qs(n) = 0,0 — m? + (1—1) + бл? 
with shape parameter бу and б). 
Lastly, 
3 3-j4j 
T(u v) = шайса) т (3.6) 
where 


Ko = бузба K, = 0j zai бізді + һубу+ааб]+ал + 2011 
К = ба eii hj9; Gris. + 26; 1151, Кз = бу1 бді. 
0.0) = @;+11(1— 07? + d —n)n + булат? 


With shape parameters Q;41 and 6,11. 
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4. THE PROPOSED ALGORITHM 


This section uses the rational cubic function and rational bi-cubic 
function to build convexity-preserving for regular surface data (2). 


4.1. Convex Curve Model 


Assume that (uj, vj) ;j = 0,1,2,..., n} represents the given convex data 
set defined over the interval [а, b] where a = ug < Uy < u; < ::: < Un = 
b and dj; j = 0,1,2, ...,n designates the derivative values at the point шу;ј = 
0,1,2, ...,n. 


Also, let 
Y Vjgi-Vj 
A;= Р 


Additionally, the following condition is fulfilled by the derivative 
values: 


A;< Ау 
а; < ҮЛЕН а; < A;< dj+1, 
о; and oj are the parameters in the interval [u;, uj +1]. 


The piecewise rational cubic function T;(u) defined in Eq. (2.1) is 
convex, if 


7, (u) > 0,uo < u < и, 
The value of T? u) is computed as follows: 


Èj- -7 08 


T”) = = (4.1) 
hj (Q; (à) 
where 
à, = —20j0jdj,, — 40;d; — Зоуоуйу + 5g;o;Á; + 40;A,, 
à = Зоі — 0/64) — 20/4; - 60)0)4) + 7оуоуАу + 20А — 
Зо/А;, 
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дз = 60/64 + 26d +1 — 3034; + 05d; + 30jÀ; – 7004; — 
à, = 300301 + 40; dj, + 20;6jd; - 5оуоуАу - 4оуду, 
as = —20;0/ Á; + 40? А, + 20j0j dj + 407 dj41. 
The denominator quantity Q;(¢) in (9) is a cube quantity and is positive 
if 
0; > 0, 0; > 0. 
In addition to this, the requisite condition for convexity in each 
subinterval is given as: 
a; = 0; j = 0,1,2,3,5. 


After numerous manipulations, we successfully procure the constraints 
on o; and o; given by 


2(d; — 4) 
oj >max10,———— (4.2) 
(-—dj,1 — 6d; + 7А,) 
2(Á; — d; 
0) > max fo, i] (4.3) 
(dj + 6d;,, — 74) 


which can be rearranged as 
2(d; — Á; 
0, — l Kj >0 
(—dj,1 — 6d; + 7А,) 
2(Aj — dj+1) | 


0, = 
(dj + 6dj,4 – 7À;) 


0} = у + max 


oy = s + max & > 0 


The above discussion is summarized as: 


Theorem 1: Assume that (uj, vj);j = 0,1,2, sig n] be a convex data set. 
Then, the convexity is sustained by the piecewise rational cubic function 
defined in (1), if in each subinterval [uj uj], j = 0,1,2,..,п the 
parameters оу, o; satisfy the conditions : 
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0; = K; + max [o = K; > 0 
К Саса ор * 
J J 


2(Á, — йы) | 
(dj + 64, — 74) 
5. CONVEX SURFACE MODEL 


tj = s max, & > 0 


Consider a convex data set (uj, v, G; = 0,1,2,...,n, 1 = 
0,1,2, suit) demarcated over the rectangular mesh (ш,ш44| X [0,051] 
fo j= 0,1,2,...,n-1,1=0,1,2,..,m—1 such that the following 
necessary conditions for convexity are settled: 


Аа ku, Gi < бала Да Na Gn S Ga 


у G; -Ĝi ~ G; -(; 
— Vrti1l jl — “j l+1 jl 
Aj m h; А A= h 
j l 


prs convexity for the surface pe defined in (3. 1) holds if all the 
boundary curves (3.3), (3. 4), (3. 5) and (3.6) are convex. 


Consequently, the boundary curve Т(и, vj) is convex if T^ (i, vj))z0ie. 


уа — )5 16, 


(2) = 
T uv) = огуз ' 


(5.1) 


where 
bo = 2020314; а + 407A; — 2010916) кы - 40716), 
Һ--20,616111-40,16%-30,6168 + 50; 10; Aj. + 401, 


b, = 30; e ЕТ Шы = 20,16 = болоб + 70j 0j + 
AE = Зоѓ Aj 


Ёз = 601616141 20516544 — 305,654 + олобу + Зоў Ау 
703191 26185 


b, = 30; 19; Ghar + 405,654, + 20510516], — 50j10j A). — 461411, 
bs = 203100, + 40А, + 205107, 6], + Асб уд. 
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Now, for positive T^ vj), the numerator Y;.9(1 — £)9^!£!b, > 0, 
If 
b, > 0,1 = 0,1,2,3,4,5 


After frequent calculations, following constraints were inferred, 


ш> max fogta) 


Gu +66", 1-74) 


Correspondingly, T (и, 10.1) is convex if T (и, 0+1) > 0 iee., 


zal vss 
h(Q4G)? ^ 


T (y, ты) = (5.2) 


where 
бу = 205 1410; д+1А +1 + 407 ае 207401143404 — 
Co = 407 1419,1414),141 Оўл+1^\]1+1 Оў1+10]1+16]+11+1 
2 fu 
407 LA 6j 141) 
Ë — ғы ғы ғы 
с = —20]1+101+1б]+11+1 = 40 j141Gj 141 Eu 3071419141 Gj 141 + 
5Qj,1419 141 уі + 405 141A; 141, 
A unm 2 хи хи хи 
C2 = 30/141 61 m 25,419,141 jai m 207.1415) 141 Е 
ЕА) - әй "M 
60 5,141 141 Gf ii + 70 p41 G41 Аі + 205 aA ipea — 30f1414j, 1+1 
A = рм рм 2 хи 
Сз = 60;1+1031+10у+11+1 + 20j 1G ja — З0у аі за + 
a "p" 2 Р 
Ojix10j I1 8j 14i + 30F1414j 1+1 703103 iau 20141, 
Жа а ғы ғы ғы 
ёд = 30j 104a jaa 4616) + 201613161 
50]1+101+1А\у1+1 - 40j AB 141, 
m $. X a. x 2 pu 
5 = —20j1+10j1+14j 1+1 + Фо] 1А +1 + 20) 1410141 G41 + 
2 Fu 
Абы біз: 
5 5-lzla 
For the numerator Y7-9(1 — £)? £t ó, = 0, 


We must have 
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ó, > 0,l = 0,1,2,3,4,5. 
which yields, 


2 Giza Ара 
0; 141 > max 40, (бише) | , 
> IG ETE, 
]+11+1 1+1 Plt. 


Äi- 
Оуу > max fo, elo 


GU +664 аа -7Aj 41) 
Identically, Т(и;, v) is convex if TPO (a v) > 01е., 


= ; 5.3 
(030002 de 


TO (uv) = 
where 
fy = 207,8; Aj. + 40710. — 20016 — 407 G1, 
f, = —20,,05,651,44 — 405165, — 305,0), 651 + 505,0; À + 40; À у, 
f, = 36), Gir = 05185165144 — 2051671 — 6616167 70510; + 
20; Àj, — 30141, 
fs = 60106), на 24167,4 — 30310, + Ojaja + 307 Aj, 
70;,0;,À;, — 20А, 
f, = 305105161144 + 495,651, 205105161 — 50б — 45; Àj, 
f; = —20;,0f Àj, + 407,Aj + 205,07, 6], + 495,651. 

Likewise, 15.9(1 — £)*7 2f; > 0, 
if 
f; > 0;j = 0,1,2,3,4,5. 


This provides the following constraints, 


C d +7Aj, 
2(Aj1 ) 
= 4-1 
Орт > тах Т Er 
А 1+1 
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Eventually, Т(и+, v) is convex if Т Gus v) = Oie., 


Yeo ЗЕ; 


Qu _ djzol 
y бы) оуу 


(5.4) 


where 
ŝo = 202.1 16: Атна 2054103231101. 

0 = 4074119741 147411 27 41,14)41,1 2741197411 5744141 — 

~2 
401.4, 1. END 

à — = = = Ат = = Av 
51 = — 2641154116); +1141 — 40уғ116)11 - 30 j411F 41145411 + 
5Q 4119 si eil + 40 xi eade 

A 332 ~ = Av Е 
52 = 302110) +11+1 — 54119} бунын = 203116. +11 
60 j4110j41, KH {+11 + 70j+110j+1, TUE, + 207411 Aj+41 — 307. aa 
А — Er = ~2 
5з = 604115) 4116); +1141 + 26)4116)41141 — 3031167, er 
брал унала + 305a 7 Tja — 26-1 Аалы» 
S4 = 30 j4115)411G/. $1141 + RETIUM t 204119 411G/. +1,1 — 

5Q 41,19 sii ei - 48; xi et 
55 = —20;44,40, НЕТ +11 + 46, oP Ay +1,1 + 2054119, 7 м6); + + 
46; Feat pea ltt 
Likewise, У (1 — Э БЕЗЕУ > 0, if 
§ > 0;j = 0,1,2,3,4,5 


This provides the following constraints, 


= 2(@1-Ау+ы) 
б}+ы > Max 40, 6—6, 
(- Geile 667,, +78 +) 
2(A jet, 1- GPa is) | 


Qj+1 > max 10, = 
(бї, [E66 1+1 -7Аҙғы) 
This discussion is briefed as: 


Theorem 2: The visualization of 3D convex data as convex surface for Bi- 
cubic partially blended rational function in (3) can be rendered if in each 
rectangular mesh I= [uj ш+] х [0,1011], the parameters 
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Qj» GF Озал, Gis Qj lo Gr Ojo бүза, assure the following 
constraints: 


0j, > max fo, 


2(0%-А,!) 2(Аы-@ ү) 
) ° 


— s 0ji > шах 40, — —— 
(- Gu 1-66 s (& +6% , 1-74 


20691-4, на) | 


0; > max 40 
j,l+1 , хи ~ ° 
| (- uaa 6G; ZEE 


Avia = G 
Qj +1 > Max fo, ae 


u pu ^ 2 
Grits +6Gj, 1,1+1 E 


= (т) = 2(Хи- буа) 
Oj > max fo, aS} Oji > max fo, CEE 


j,l+1 
2(67,4 )-A; 
= +11 TELL 
Oj+1,1 > Max Que) | , 
(761a 66% 7 Ағы) 


2(Á ssi б, 1) | 


Qj+1 > max 10, ; 
"(ау [+1 159687, 1 7Ë au) 


Adjusting the above constraints as: 


2j — Аи) |, 
о; = a; + max ——v M, —.— a; > 0 
RU EIU à 
2 GA 
Oji = bj, + max fo, МОУ Аш) | bj, >0 
"(ба баў ы) 
(С — A; 
Oj 41 = Сул + max [o em йы ај 0 Cji > 0 
| j+1,1+1 6G, 4 + 7А) 
2(Aj ria — ба) | 
š = d; , + max ТЕП m mw. d, > 0 
на d | (бу 1+1 + 6 ааа 75,144) ü 
2(6j, — | 
gGim04g-q4magx10,2——————— —-— ej; > 0 
d и Ë = j, Жат ja +74 м) d 
2(Kj, 27 E | 
0;; = gj, + max Ша eT к ТТТ |) ы >0 
Qj = Bj fo, TE ӨС _ 7K j Jj, 
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2(G7 1,⁄ Қы 
Oj +11 = mj, + max 0, ae | mij, > 0 
(Ga I+1 77 607411 + 7Aj+11) 
2(41- ) 
= Ga 1+1 
бұғы = Yt max, ая пуд > 0 
(G j+1,1 + 66 te Т 7Kj +11) 


6. RESULTS AND DISCUSSIONS 


This section describes algorithm for the proposed shape preserving 
techniques in section 4 and verify its effectiveness through graphical 
results in the subsequent numerical examples. 


6.1. Algorithm. Convex Curve Model 
Step 1: Compute (n + 1) convex data points (uj, vj);j = 0,1,2, -— 
Step 2: Approximate the derivatives A; at the nots using Arithmetic 
mean method. 
Step 3: Using Theorem 1, evaluate free parameters o; and о}. 
Step 4: Insert the above values in Eq. (4.1). 


Example 1. The data in Table 1 encloses convex data set. MATLAB 
programming software packet was used to implement constraints on shape 
parameter in order to visualize convex curve. Numerical results were 
recorded in Table 5. Figure 1 clearly reflects the shape characteristic of 
convexity. 


Table 1. A Convex 2D Data Set 


Uj -3 —2 —1 0 1 2 3 4 5 


Vi 17 10 5 1.8 1 1.8 5 10 17 


Example 2. Another convex data set was considered to implement the 
proposed algorithm using MATLAB software packet and results are 
shown in Table 6, while Figure 2 clearly depicts convexity of data. 


Table 2. A Convex 2D Data Set 


Uj 1 2 3 4 4.9 6 
Vi 4.11 9.39 20.3 54.2 163.89 4381 


Scientific Inquiry and Review 


х УК 


Volume 7 Issue 2, 2023 


Ibraheem et al. 


v-axis 


v-axis 
Figure 1. Convexity Preserving Curve of Convex 2D Data Set in Table 1 
Using Convex Rational Cubic Spline Function 


250 


v-axis 


1 1.5 2 25 3 35 4 45 5 55 6 
u-axis 

Figure 2. A Convex Curve of Data Set taken in Table 2 Using Convex 

Rational Cubic Spline Function 


Algorithm 2. Convex Surface Model 


Step 1: Compute (n + 1) x (m + 1) convex data point 
points ((uj,vj, G; i): = 0,1,2,...,n; L= 0,112, ., m }. 
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Step 2: Approximate the derivatives 


Arithmetic mean method. 


G?, and Gat the nots using 


J 


Step 3: Using Theorem 2, evaluate free parameters Oji, jl, Qj141; 


бізі» Oj, Oj1- Qj+11 and бууу. 


Step 4: Insert the above values in (3). 


Example 3. A convex 3D data set is prescribed in Table 3 which is created 
by the function particularized in the following Eq. (6.1). Convexity is 
conserved by utilizing convex bi-cubic rational function with the graphical 
representation illustrated in Figure 3. Numerical results are displayed in 
Table 7. The uz-view and the vz-view of the convex surface is given in the 
Figure 4 and Figure 5, respectively. 


G(u,v) = u? + v? + 0.1. (6.1) 

Table 3. A Convex 3D Data set 

uv -5 4а 3 2 4 0 1 2 3 4 5 6 
5 501 411 341 291 261 251 261 291 341 41.1 501 611 
4 41.1 321 251 201 Пі 161 171 201 251 321 411 52.1 
3 341 251 181 131 101 91 101 131 181 251 341 45.1 
2 291 91 131 81 51 41 51 81 13.1 201 291 401 
4 261 171 101 51 21 11 21 51 101 171 261 371 
0 251 161 9.1 41 11 01 11 41 91 161 251 361 
1 261 171 101 51 21 11 21 51 101 171 261 371 
2 291 201 131 81 51 41 51 81 131 201 291 401 
3 341 251 181 131 101 291 101 131 181 251 341 4541 
а 411 321 251 201 171 161 171 201 251 321 411 521 
5 501 50.1 41.1 341 291 261 251 261 291 341 41.1 501 
6 611 521 451 401 371 361 371 401 451 521 611 721 
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Figure 3. A Convex Surface produced from the Convex Data Set in Table 
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Figure 4. A uz-view of the Convex Surface 
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80 


z-axis 


Figure 5. A vz-view of the Convex Surface 


Example 4. Another convex 3D data set in Table 4 is made up by the 
function given by the subsequent equation (17). The convex rational bi- 
cubic functions help to preserve the convexity of the data which can be 
visualized in Figure 6. Figures 7, and Figure 8 give the uz-view and vz-view 
of the convex surface, respectively. Numerical results are shown in Table 
8. 


(и, v) = (u? + v? — 1)? + 0.2. (6.2) 

Table 4. A Convex 3D Data set 

u/v 0 1 2 3 4 5 6 7 
0 1.2 0.2 9.2 64.2 225.2 576.2 1225.2 2304.2 
1 0.2 1.2 16.2 81.2 256.2 625.2 1296.2 2401.2 
2 9.2 16.2 49.2 144.2 3612 784.2 1521.2 2704.2 
3 64.2 81.2 1442 289.2 5762 1089.2 1936.2 32492 
4 225.2 2562 3612 576.2 9612 1600.2 2601.2 4096.2 
5 5762 6252 7842 1089.2 1600.2 2401.2 3600.2 5329.2 
6 12252 1296.2 1521.2 1936.2 2601.2 3600.2 5041.2 7056.2 
7 2304.2 2401.2 2704.2 3249.2 4096.2 5329.2 7056.2 9409.2 
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Figure 7. A uz-view of the Convex Surface 


UMT— 21 


School of Science жа 
© 
Volume 7 Issue 2, 2023 Pun 


Shape-Preserving Curve... 


z-axis 


о 1 > 3 et 5 6 = 8 
v-axis 


Figure 8. A vz-view of the Convex Surface 
7. NUMERICAL RESULTS 


This section presents the numerical results when proposed algorithm was 
implemented on 2D convex data set in Table 1 and Table 2. Similarly, for 
3D convex data in Table 3 and Table 4, numerical results are given in Table 
7 and Table 8. MATLAB software packets were used for programming and 
visualizing convex curves and surfaces. 


Table 5. Numerical Results Parallel to Figure 1 


i 1 2 3 4 5 6 7 8 9 
d, -8.000 -6.000 -4.100 -2.000 0 2.000 4.100 6.000 8.000 
0; 1 1 1 1 1 1 1 1 
Oi 1 1 1 1 1 1 1 1 

Table 6. Numerical Results for Figure 2 
i 1 2 3 4 5 6 
di 2.4650 8.0950 22.4050 77.8889 185.5798 319.3540 
0; 1 1 1 1 1 
Oi 1 1 1 1 1 
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Table 7. Numerical Results for Figure 3 


u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 


Numerical results for би 


-5 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 -10.0 
-4 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 -8.00 
-3 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 -6.00 
-2 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 -4.00 
-1 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 -2.00 
0 0 0 0 0 0 0 0 0 0 0 0 0 

1 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 
2 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 
3 6.00 6.00 6.00 6.00 6.00 6.00 6.00 6.00 6.00 6.00 6.00 6.00 
4 8.00 8.00 8.00 8.00 8.00 8.00 8.00 8.00 8.00 8.00 8.00 8.00 
5 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 10.0 
6 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 12.0 


Numerical results for б}, 


-5 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
-4 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
-3 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
-2 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
-1 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
0 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
1 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
2 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
3 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
4 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
5 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 
6 -10.0 -8.00 -6.00 -4.00 -2.00 0 2.00 4.00 6.00 8.00 10.0 12.0 

Numerical results for Q; , 

-5 1.00 1.00 1.00 1.00 1.000 1.000 1.000 1.000 1.000 1.000 1.000 AES 
-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mm 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
6 = a 

Numerical results for o; , 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
d 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 с 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 === 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mm 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
6 


Numerical results for Qj 141 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 

1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 с 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 === 
6 == == 

Numerical results for G; 141 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2. 
-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EM" 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
6 

Numerical results for о; 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 === 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mm 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
6 TN — 

Numerical results for 9; , 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 == 
-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 с 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 === 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mm 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
6 

Numerical results for Qj 41, 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 zu 

-4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v -5 -4 -3 -2 -1 0 1 2 3 4 5 6 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 

1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 с 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
6 a. — 

Numerical results for 07411 

-5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 "EN 
d 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mm 
-3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
-1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 __ 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v -5 -4 -3 -2 1 0 1 2 3 4 5 6 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ME 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 NES 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 X 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
6 

Table 8. Numerical Results for Figure 4 
u/v 0 1 2 3 4 5 6 7 

Numerical Results for Gul 

0 -0.0060 -0.0060 -0.0060 -0.0060 -0.0060 -0.0060 -0.0060 -0.0060 

1 0.0040 0.0080 0.0200 0.0400 0.0680 0.1040 0.1480 0.2000 

2 0.0320 0.0400 0.0640 0.1040 0.1600 0.2320 0.3200 0.4240 

3 0.1080 0.1200 0.1560 0.2160 0.3000 0.4080 0.5400 0.6960 

4 0.2560 0.2720 0.3200 0.4000 0.5120 0.6560 0.8320 1.0400 

5 0.5000 0.5200 0.5800 0.6800 0.8200 1.000 1.2200 1.4800 

6 0.8640 0.8880 0.9600 1.0800 1.2480 1.4640 1.7280 2.0400 

7 1.2940 1.3220 1.4060 1.5640 1.7420 1.9940 2.3020 2.6660 
Numerical Results for б? и 

0 -0.0060 0.0040 0.0320 0.1080 0.2560 0.5000 0.8640 1.2940 

1 -0.0060 0.0080 0.0400 0.1200 0.2720 0.5200 0.8880 1.3220 

2 -0.0060 0.0200 0.0640 0.1560 0.3200 0.5800 0.9600 1.4060 

3 -0.0060 0.0400 0.1040 0.2160 0.4000 0.6800 1.0800 1.5460 
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u/v 0 1 2 3 4 5 6 7 
4 -0.0060 0.0680 0.1600 0.3000 0.5120 0.8200 1.2480 1.7420 
5 -0.0060 0.1040 0.2320 0.4080 0.6560 1.000 1.4640 1.9940 
6 -0.0060 0.1480 0.3200 0.5400 0.8320 1.2200 1.7280 2.3020 
7 -0.0060 0.2000 0.4240 0.6960 1.0400 1.4800 2.0400 2.6660 
Numerical results for о; 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _— 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 AN 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2. 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 I 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
7 — --- == >= === — — a 
Numerical results for 0; 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 en 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 T 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 "M 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ___ 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
7 === — — — = == — c 
Numerical results for 0/11 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ЕЕ 
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u/v 0 1 2 3 4 5 6 7 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 M 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ___ 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _— 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _— 
7 — — — === — — — = 

Numerical results for o; 141 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2. 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 == 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2. 
7 == — — — — === — шіт 
Numerical results for 0) | 

0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mE 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 MEN 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EM 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ES 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 5. 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 
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u/v 0 1 2 3 4 5 6 7 
7 --- ---- — — — —. —— ЕРЕ, 
Numerical results for 0; 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _ 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 E 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2 2 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 "n 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 БЕН 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ERN 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _ 
7 — — E === — — — === 
Numerical results for @j41,1 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ES 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 _ 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 = 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 mn 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 a 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 NM 
7 ——— —— — === == — — —— 
Numerical results for 0411 
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ___ 
1 1.000 1.000 1.000 1.000 1.000 1.000 1.000 MET 
2 1.000 1.000 1.000 1.000 1.000 1.000 1.000 un 
3 1.000 1.000 1.000 1.000 1.000 1.000 1.000 ии 
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u/v 0 1 2 3 4 5 6 7 
4 1.000 1.000 1.000 1.000 1.000 1.000 1.000 EN 
5 1.000 1.000 1.000 1.000 1.000 1.000 1.000 — 
6 1.000 1.000 1.000 1.000 1.000 1.000 1.000 2 
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8. CONCLUSION 


In the current research, shape preserving rational schemes were devised 
for convex curve and surface data. A C1 rational cubic piecewise function, 
composed of cubic numerator and quadratic denominator with two free 
parameters o; and 0), was used to retain the intrinsic shape characteristic of 
2D convex data. Mathematical derivation was performed to derive the data 
dependent constraints on free parameter and results were summarized in 
Theorem 1. MATLAB software packet was used to implement the 
developed scheme and visualize 2D convex data. Furthermore, the function 
was then advanced to partially blended rational bi-cubic function with eight 
free parameters for the 3D convex data positioned over rectangular grid. 
These eight parameters оу, 051, оза, 05115 9а, OQj1 G41 and Oj41,1 
by means of mathematical derivation, were conditioned to preserve the 
inherent curvature of 3D data and findings were summarized in Theorem 2. 
The proposed algorithm controls the shape of visual model locally. Convex 
surfaces were realized with the help of MATLAB software packet. The 
suggested shape-preserving techniques in this study are applicable to both 
data and derivatives and do not alter basis functions. The proposed scheme, 
however, does not use free parameters which provides the users with 
minimal flexibility to change the shape of the data as they desire. Work on 
curvature preserving scheme for scatter data is in process by the authors. 
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ABSTRACT 


The current study examined the toxicity of Roundup, which is a herbicide 
containing glyphosate that is commonly used in agriculture, on 
Oreochromis niloticus, a type of freshwater fish. The researchers 
investigated the effectiveness of using antioxidant enzymes as biochemical 
biomarkers to assess the potentially toxic impact of Roundup on the fish. 
The sub-lethal concentration of Roundup in tilapia is 1.05 mg/l for 96 hours 
of exposure. The fish were exposed to sub-lethal concentration (36.8 to 
46.9 mg/L) of Roundup for a period of 24, 48, 72, and 96 hours, and their 
antioxidant enzyme activities were measured, respectively. According to 
the findings, Roundup exposure led to a noticeable rise in the levels of 
superoxide dismutase, catalase, and peroxidase activities in both the liver 
and gill tissues of the fish, which suggested the occurrence of oxidative 
stress. These changes in antioxidant enzyme activities were concentration- 
dependent, suggesting that Roundup exposure can induce oxidative damage 
in fish. The study concluded that Roundup affects the enzymatic activity of 
both organs in a concentration- and time-dependent manner. Catalase and 
peroxidase activities were generally inhibited, while superoxide dismutase 
activity was affected differently depending on the organ, concentration, and 
time of exposure. 


Keywords: biochemical, biomarkers, enzyme, fish, gills, liver 
1. INTRODUCTION 


Due to the fact that approximately 7096 of Pakistan's population relies 
on the agricultural sector, there has been a significant increase in the use of 
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pesticides over the last 30 years to meet the growing demands of the 
population [1]. The accessibility and affordability of these pesticides have 
led to their indiscriminate use, resulting in significant ecological changes in 
major freshwater resources such as water, sediment, and fish [2]. 
Glyphosate, an organophosphorus herbicide, is currently the most 
commonly used herbicide on a global scale [3]. It is marketed under 
numerous brands and commercialized in over 100 countries worldwide [4]. 
The primary toxic component of formulated glyphosate-based products for 
aquatic organisms are surfactants like POEA found in Roundup. Several 
previous studies have suggested that pesticides can cause alterations in the 
water present in the ponds surrounding the cultivation regions, as well as 
rainwater and groundwater [5] through mechanisms such as volatilization, 
spray drift, surface runoff, chemical-biological degradation, leaching, and 
sorption [6]. 


Antioxidant enzymes, including CAT, POD, and SOD, are important in 
fish because they protect the fish from oxidative stress caused by factors 
such as pollution and diseases. These enzymes help to neutralize harmful 
free radicals and reactive oxygen species, which can damage cells and 
tissues. Additionally, these enzymes play a crucial role in the overall health 
and survival of fish, as oxidative stress can lead to a range of negative 
effects, including reduced growth rates, decreased immune function, and 
increased mortality. Therefore, understanding the role of antioxidant 
enzymes in fish can help to inform efforts to protect and maintain healthy 
fish populations [7]. 


There is limited information available in the previous literature 
regarding the harmful effects of glyphosate on organisms exposed to it in 
their natural environment. However, certain synthetic agrochemical 
compounds like glyphosate may become more toxic to non-target species 
due to differences in their bioactive properties. These compounds often 
cause sub-lethal damage to tissues, disrupt endocrine function, induce 
immunotoxicity and neurotoxicity, and alter the metabolism of the affected 
animals [8]. According to Naz et al. [9], industrial and agricultural water 
sources release a range of toxins into both rivers and canals, which then 
accumulate in sand. This accumulation of toxins in aquatic environments 
has emerged as a global issue, leading to health issues for both humans and 
aquatic animals [10]. Nile tilapia (Oreochromis niloticus) is a popular 
commercial fish species and is often used as a bio indicator to assess the 
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quality of aquatic ecosystems, particularly for pollutants such as herbicides 
[11]. It is also considered a suitable biological model for toxicological 
studies due to its sensitivity to diazinon, an organophosphate pesticide 
herbicides [12]. Aquatic organisms exposed to pesticides often experience 
oxidative stress, which is caused by the production of ROS in response to 
toxicants [13]. The effects of lufenuron exposure on Oreochromis niloticus 
were examined through histopathological analysis of various organs, 
including the gills, kidneys, livers, brain, and hearts. The findings revealed 
severe alterations and indicated that lufenuron caused DNA damage and 
oxidative stress in the fish. The study aimed to investigate the mechanisms 
behind lufenuron's toxicity, as well as to raise concerns about the use and 
disposal of lufenuron in terms of its potential to cause environmental 
pollution. The study also proposed measures to mitigate contamination in 
aquatic populations [14]. The comet assay, also known as Single Cell Gel 
Electrophoresis (SCGE), as a genotoxicity technique that detects DNA 
damage caused by alkylating chemicals, intercalating agents, and oxidizing 
agents. The comet assay is a useful method for monitoring genotoxicity in 
aquatic environments, where fish are often used as test organisms to detect 
DNA damage resulting from the ingestion of mutagens and pro-mutagens 
in freshwater and saltwater. This technique has also been employed to 
determine the genotoxic potential of water resources such as rivers and 
lakes. The purpose of the current study was to investigate the adverse effects 
of lufenuron exposure on various health biomarkers of Oreochromis 
niloticus [14]. 


The objective of the current study was to examine the Roundup 
exposure affects oxidative stress in fish by assessing variations in the 
activities of antioxidant enzymes. The aim wato obtain an understanding of 
the possible adverse impacts of Roundup exposure on the well-being of 
freshwater fish and the general health of aquatic ecosystems. The aim of the 
research is to increase our knowledge about the dangers related to the usage 
of herbicides that contain glyphosate in freshwater ecosystems. 
Additionally, it seeks to lay the groundwork for more extensive 
investigations into the lasting consequences of Roundup exposure on 
aquatic creatures such as fish. 


2. MATERIAL AND METHODOLOGY 


At the laboratory of Government Sadiq College Women University in 
Bahawalpur, Pakistan, an experimental study was conducted, adhering to 
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the ethical principles for laboratory animal usage set by the Ethical 
Committee of Government Sadiq College Women University Bahawalpur. 


2.1. Sample Collection and Management 


Thirty freshwater fish of the species Oreochromis niloticus, with a 
weight range of 75-80g, were procured from the nearby Fish Seed Hatchery 
in Bahawalpur. The fish were fed with the commercial feed palettes. These 
fish were transferred to the Zoology department testing laboratory of 
Government Sadiq College Women University Bahawalpur in plastic bags 
equipped with enough oxygen. The fish were subsequently relocated to 
glass tanks filled with 100 liters of water. The experiment employed five 
100L tanks for each concentration, with each tank measuring 36 inches in 
length and 12 inches in width, housing six fish in each tank. Prior to 
commencing the experiment, the fish were given a period of 15 days to 
acclimate. To ensure that the fish had enough oxygen, continuous air pumps 
were utilized. The temperature of the aquarium water was kept constant, 
and measurements of the water's physical properties were taken before and 
after the study, as outlined in Table 1. The study was carried out in 
compliance with the ethical standards set forth by the Ethical Committee of 
Government Sadiq College Women University Bahawalpur. 


Table 1. Physical Parameters of Water 


Parameters Values 

Dissolved oxygen (mg/L) 6.70+0.1 
Electrical conductivity at 

25*(umhos/cm) 393.10--2.2 
РН 6.89+0.05 
Total dissolved solid (mg/L) 153.89--3.9 
Total hardness (CaCO3, mg/L) 174.59+2.1 
Water temperature (°C) 23.80+0.3 


2.2. Experimental Design 


During the acclimation period, 30 fish [14] were distributed randomly 
into the control group (without exposure) and four experimental groups, 
which were subjected to different concentrations of Roundup original (46.9, 
44.4, 40.0, and 36.8 mg L- 1). The selection of these concentrations was 
informed by previous research, which examined the impact of acute 
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exposure to Roundup® at comparable concentrations in other fish species 


[15]. 
2.3. Evaluation of Antioxidant Enzymes 


After exposure, fish from each group were euthanized at different time 
points (24, 48, 72, and 96 hours), and their gills and liver were extracted. 
Briefly, for removal of erythrocytes, all the organs were washed using 
phosphate buffer solution having 6.5 PH. After that all the organs were 
triturated using 0.2 mL normal saline and then homogenate was prepared 
(1:4 w/v) by adding cold buffer. The homogenate from each tissue was 
centrifuged (10000 rpm) at 40C for 15 minutes. After centrifugation, the 
supernatants were removed and placed at -800C for enzyme profiling. By 
using UV-visible spectrophotometer, the amount of (reactive oxygen 
species) ROS in the gills and liver of each fish was measured by different 
antioxidant enzymes, which were such as peroxidase [16], catalase [17], and 
superoxide dismutase (SOD) were measured according to previous 
protocols [16]. Additionally, the gill tissue homogenates from individual 
fish, which were used to quantify the level of antioxidant enzymes, such as 
catalase, superoxide dismutase (SOD), and peroxidase described by [16, 
17). 


At the Zoology Lab of the Department of Zoology, Government Sadiq 
College Women University, Bahawalpur, 36 fish samples were carefully 
dissected for the laddering assay. The dissection process was performed 
with great care to avoid contamination and the fish muscles were preserved 
in plastic bottles containing 0.9% saline solution. The bottles were properly 
labeled and stored at a temperature of -20°C, with plans to process the 
samples within two weeks. It is advisable to store muscle tissues at a 
temperature of -20°C or lower for future use, as this can enhance their 
usability [18]. 


2.4. DNA Extraction and Laddering assay 


The Phenol Chloroform Isoamyl method was used to extract DNA from 
muscle tissue samples weighing 100-200 milligrams [18]. The samples 
were homogenized and mixed with 1 ml of lysis buffer containing 1 M Tris, 
5 M NaCl, 0.5 M EDTA, 10% SDS, and 20 ul of Proteinase K. The mixture 
was then incubated overnight in a water bath at 50 degrees Celsius. On the 
following day, 1 ml of PCI (phenol chloroform isoamyl alcohol) was added 
in the prepared solution and the contents were centrifuged at 15,000 rpm for 
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10 minutes. The supernatant was transferred to fresh tubes and 20 ul of 
glycogen blue and 1 ml of chilled isopropyl alcohol were added in the 
mixture. The mixture was frozen overnight at 30 degrees Celsius and 
centrifuged at 15,000 rpm for 10 minutes. The resulting DNA pellet was 
washed three times with 70% ethanol and later dried for 15 minutes at the 
room temperature. Finally, 100 ul of injection water was added to the dried 
DNA pellet and the extracted DNA was stored at -8 degrees Celsius. To 
identify any genomic DNA damage, the extracted DNA was separated at 
1% agarose gels and visualized with ethidium bromide (6 ug/ml) under 300 
nm UV light. 


2.5. Statistical analysis 


SPSS software (version 15) was utilized for statistical analysis and 
ANOVA was employed to evaluate the significance of differences between 
the Roundup-treated and untreated fish groups. The results were expressed 
as mean + SE, and a significance level of P<0.05 was adopted. 


3. RESULTS 
3.1. Antioxidant enzyme activities 


The study investigated the impact of various pesticide doses on the 
antioxidant enzymes in different organs of fish, such as the liver and gills. 
The results revealed that the levels of catalase [17] and superoxide 
dismutase (SOD) were significantly lower in the heart, gills, liver, and brain 
tissues of fish exposed to different concentrations of insecticides as 
compared to untreated fish. Conversely, fish, which were exposed to the 
same concentrations of insecticides had higher levels of peroxidase [16] in 
their gills, liver, and brain tissues than those, which were untreated. 
According to the results of the study, there was a significant decrease in the 
levels of antioxidant enzymes (POD, CAT, and SOD) in the gills of fish 
exposed to Roundup GBH at concentrations of 0.80, 1.30, 1.90, 2.50 ug/L 
mg 17! for 24, 48, 72, and 96 hours as compared to the control group. The 
fish in group-E exhibited a substantial reduction in SOD and CAT levels in 
their liver tissue after 24 hours of exposure. Meanwhile, groups D-E 
demonstrated lower SOD and CAT levels in their liver tissue after 48 and 
72 hours of exposure as compared to the control group. Furthermore, the 
concentration of peroxidase in the liver, kidneys, and gills tissues of treated 
fish in groups D-E was significantly lower (p « 0.05) than in the non-treated 
fish. All results were presented as mean + SE, and statistical significance 
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was assessed using one-way ANOVA analysis with a significance threshold 
of p<0.05. (Table 2, 3). 


Table 2. Different Antioxidant Enzymes in Liver of Fish Oreochromis 
niloticus 


Groups 
аа A (control) B(0.80ppm) C (1.30 ppm) р Е (2.50) 
(1.90ppm) 
Catalase (units/min) 
24 4.92 +0.02 6.67 + 0.02. 4.72+0.02 6.54-0.01 5.67+0.06* 
48 4.83 + 0.03 5.64 + 0.03 4.46+0.04 5.61-0.08% 5.384 0.07* 
72 4.5640.06* 5.78 +0.01 647+-003 4.91+0.04 4.62+0.06* 
96 5.61+0.05* 5.64+0.04 5.78+0.06* 4.56+0.07* — 4.902 0.03 
Superoxide dismutase (units/mg) 
24 12.98+0.11 12.9140.01 12.86+0.01 11.67+0.05 9.11+0.03* 
48 13.04+0.03 12.75+0.06 12.50+0.02 9.20+0.49* 9.00+0.05* 
72 13.07 +0.05 12.80+0.04 12.29+0.05 8.22+0.29* 8.37+0.01* 
96 11.91+0.07* 13.45+0.02 11.77+0.05 9.4140.07* 9.45+0.01 
Peroxidase (units/min) 
24 5.21+0.01 6.11+0.02 5.07+0.02 4.79+0.02* 4.23+0.01* 
48 5.20+0.05 612+0.02 2.964002 2.61+0.03* 4.11 £ 0.01* 
72 5.2340.02 6.09+0.08  2.8540.04 2.2440.04* 4.02+0.01* 
96 3.9440.02 6.07+0.05 5.1940.02 5.71+0.01 3.10+0.03 
Table 3. Enzymes in Gills of Fish Oreochromis niloticus 
Parameter/hour Groups 
4. A (control) B (0.80) pn pun E (2.50) 
24 3.23+0.00 3.254001 3.18+0.03 3.68+0.01 3.34+0.04* 
48 3425000 3.114002 317-002 2.8340.06* 3.27+0.03* 
72 3.06+0.00 2.87+0.01 247-006 3.23+40.04* 3401-0401 
96 3115000 3.364002 3.12+0.02 345-001 2.26+0.01 
Superoxide dismutase (units/mg protein) 
24 5.85+£0.04 5.73+0.04  5.63:0.0 4.89+010 5.64-0.04% 
48 579+003 5.56+0.01 5.88+0.03 4.16+0.03* 5.29+0.01* 
72 5.93+0.04 5.94+0.02 626-001 4.99-0.04% 5.48+0.11* 
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Parameter/hour Groups 

caca A (control) B (0.80) 2. pes E (2.50) 
96 5.97 + 0.03 6.45 + 0.02 5.35+0.02 5.02+0.05* 4.96+0.09* 
Peroxidase (units/min) 

24 1.42+0.01 1.29+0.01 1.22+0.01 1.314+0.04* 0134+0.02 
48 1.34 + 0.01 1.31 +0.03 1.28+0.02 1.29 € 0.05* 1.324 0.01 
72 1.37 x 0.01 1.28 + 0.02 1.26+0.01 1.25 +0.03* 1.29+ 0.02 
96 1.45+ 0.03 1.31 0.01 1.31+ 0.02 1.34+ 0.01 1.41- 0.01 


3.2. DNA Quantification 


Nano Drop technology measures the absorbance of the light by a DNA 
sample to determine its concentration and purity. A quick and easy method, 
which requires only a small amount of the sample. The technique is 
commonly used in research and diagnostic laboratories for DNA 
quantification [19]. The Nanodrop technique was employed to measure both 
the quantity and quality of genomic DNA using a Spectrophotometer 
(Agilent 8453, California, USA). The absorbance values of the DNA 
sample, measured at 260/280 nm, were recorded for 5 seconds (Figure 1) to 


determine the quantity and purity of the sample. 
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Figure 1. Absorbance Peaks Showing Good Quality of DNA 
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The quality and quantity of genomic DNA were assessed using the 
Nanodrop technique and a spectrophotometer (Agilent 8453, California, 
USA). The absorbance values at 260-280 nm wavelength were used to 
evaluate the DNA quality and DNA concentration values (ng/ul) and the 
ratio of 260/280, which were noted for purification purposes. The following 
equation was utilized to calculate the DNA concentration: volume (g/ml) = 
OD 260 x 100 (dilution factor) x 50g/ml/1000. The results were acceptable 
only if the ratio of OD 260/OD 280 was 1.8-2.0, indicating good DNA 
quality. Table 4 summarizes the obtained data. 


Table 4. DNA Quantification using Spectrophotometer (Agilent 8453, 
California, USA) DNA damage 


24 days Insecticide concentration 260/280 ratio DNA concentration 

Gills 1.49 44.9ug/ml 
0.80u g/L 

Liver 1.28 44.65ug/ml 

Gills 1.28 4lug/ml 
1.30u g/ml 

Liver 1.95 24.1u g/ml 

Gills 1.67 2425ug/ml 
1.90u g/ml 

Liver 1.65 23.45ug/ml 

Gills 1.67 2425ug/ml 
2.50g/ml 

Liver 1.65 23.45u g/ml 

48 days Insecticide concentration 260/280 ratio DNA concentration 

Gills 1.58 22.7ug/ml 
0.80u g/L 

Liver 1.47 2lug/ml 

Gills 1.41 20.5ug/ml 
1.30u g/ml 

Liver 1.71 18.75ug/ml 

Gills 1.24 17.65ug/ml 
1.90u g/ml 

Liver 1.13 16.35ug/ml 

Gills 1.41 20.5u g/ml 
2.50g/ml 

Liver 1.58 22.7ug/ml 

72 days Insecticide concentration 260/280 ratio DNA concentration 

Gills 0.80u g/L 1.54 16.25ug/ml 
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72 days Insecticide concentration 260/280 ratio DNA concentration 
Liver 1.47 14.65ug/ml 
Gills 1.60 14.35u g/ml 
LACE 1.30u g/ml 
Liver 1.40 133p g/ml 
Gills 1.37 13ug/ml 
1.90u g/ml 
Liver 1.26 12.3ug/ml 
Gill 1.44 16.26u g/ml 
2.50g/ml 
Liver 1.46 14.66u g/ml 
96 days Insecticide concentration 260/280 ratio DNA concentration 
1.54 16.24ug/ml 
Gill 1.46 14.63u g/ml 
— 0.80ug/L 
Liver 1.68 14.34u g/ml 
Gill 1.47 13.3ug/ml 
== _, 1.30u g/ml 
Liver 1.36 13.00u g/ml 
Gill 1.25 12.10u g/ml 
rud 1.90u g/ml 
Liver 1.45 16.21ug/ml 
Gill 1.47 14.61ug/ml 
== 2.50g/ml 
Liver 1.46 14.66u g/ml 
4. DISCUSSION 


Glyphosate, a primary component in numerous herbicides, is 
extensively employed in agriculture industry to eliminate unwanted weeds 
and in non-agricultural areas for vegetation clearance. Its usage can lead to 
contamination of nearby water bodies through irrigated food crops, which 
can cause various abnormalities in the organisms inhabiting the polluted 
water. Due to its frequent use and specificity, glyphosate has been found to 
induce several toxicological problems in exposed species, including 
humans. Organisms' antioxidant systems play a crucial role in responding 
to oxidative stress by generating antioxidant metabolites and enzymes. 
Catalase is a significant enzyme among these enzymes, as it helps in 
preventing the harmful impacts of ROS on cells [20]. In the current 
investigation, the liver tissue levels of lipid peroxidation, which is the 
primary immune organ, exhibited a significant decrease (P < 0.05) in groups 
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exposed to 220 mg/L of GBH. Moreover, the catalase activities in liver 
tissues also increased. Previous research has reported significant alterations 
in oxidative stress in fish upon exposure to GBH [21]. The current study 
found that fish exposed to glyphosate had increased levels of ROS and 
TBARS, while their levels of GSH and total proteins decreased, 
respectively. The liver tissues of treated fish showed these effects 
throughout the experiment. Additionally, the kidneys of fish in groups D-E 
showed significantly higher levels of oxidative stress and ROS contents 
after the 4th day of exposure as compared to unexposed fish [22]. The levels 
of lipid peroxidation in liver tissues of fish exposed to GBH exhibited 
diverse outcomes, either increasing, decreasing, or remaining stable, 
depending mainly on the fish weight (or age) and exposure time, similar to 
the levels of catalase. Such oxidative conditions have been observed in 
various fish species upon exposure to different pollutants, including GBH, 
as per earlier studies. In this study, the DNA fragments displayed irregular 
and stepped smears due to the inaccurate process of DNA cleavage that 
occurs during apoptosis. This result is consistent with earlier findings [23] 
conjecture that more than one mechanism could be responsible for the 
fragmentation patterns leading to cell death and the cumulative effect of 
these events can trigger apoptosis and necrotic cell death. The ends were 
employed to assess the intricate genotoxic effects of the insecticide at the 
DNA level of many fish species under study. Various tissues are commonly 
used to determine DNA damage, namely , gills and liver [24]. In this study, 
the liver -detoxification organ- and the gills were selected as organs that are 
directly exposed to environmental contaminants. Gills and liver, being 
metabolically active organs that have higher metal accumulation than other 
organs, share this characteristic. Metal ion complexes formation can cause 
DNA damage. During apoptosis, DNA cleavage occurs at sites between 
nucleosomes, leading to DNA damage, which is a typical feature [25]. 
Moreover, chromatin is broken down by proteases and endonucleases into 
a smearing pattern, which is an indication of necrosis when proteases break 
down the histones and expose the entire length of the nucleases. 


4.1. Conclusion 


By analyzing the activities of antioxidant enzymes as biochemical 
biomarkers in freshwater fish Oreochromis niloticus (glyphosate-based 
herbicide) exposed to Roundup, it was concluded that the exposure to 
Roundup resulted in oxidative stress as well as DNA damage in the fish. 
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Therefore, the current conclusion was supported by the changes noticed in 
the activities of antioxidant enzymes. The study found a significant increase 
in oxidative stress by the activities of superoxide dismutase (SOD), 
catalasend peroxidase. The active ingredient in Roundup and glyphosate, 
can generate reactive oxygen species, which can damage the DNA 
molecules. This damage can lead to mutations and genetic instability, 
potentially increasing the risk of cancer and other diseases. Moreover, this 
sort of damage can lead to reduced growth, reproduction, and survival of 
fish populations. Additionally, glyphosate can also affect the aquatic food 
chain, potentially affecting the health of other species, which rely on fish 
as a prime food source. These findings emphasized the need for proper 
regulation and management of herbicide usage in aquatic ecosystems to 
minimize the impact on fish and other aquatic organisms. Thus, it is crucial 
to monitor the use of glyphosate-based herbicides іп freshwater 
environments and take measures to reduce their potential negative impacts 
on aquatic organisms. Exposure to Roundup herbicide led to changes in the 
DNA concentration and 260/280 ratio in fish organs such as the liver and 
gills. The gills showed higher vulnerability to Roundup herbicide exposure 
than the liver. This suggested that Roundup herbicide exposure could cause 
DNA damage in fish organs. Hence, further research is required to evaluate 
the long-term effects of Roundup herbicide exposure on fish populations 
and the broader ecosystem. For this purpose, measures should be taken to 
reduce Roundup herbicide usage and minimize its potential negative 
impacts on aquatic life. Further research can also comprehensively 
understand the long-term effects of Roundup exposure on fish and other 
aquatic organisms. 
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ABSTRACT 


Human beings heavily rely on the sense of sight, as sight mainly controls 
approximately 8096 of our perception. During the COVID-19 pandemic, 
students had to take their online classes by using different electronic 
devices. By not following the proper instructions to use these devices, 
primary students faced eyesight issues during the COVID-19 pandemic. 
Low Vision is a commonly diagnosed issue in ophthalmology in Pakistan. 
Therefore, the current study is a survey-based research. For this purpose, a 
population of 100 primary students was gathered who took online classes 
during the pandemic and developed eyesight issues after their online 
classes. Based on statistical analysis, the logistic regression model was used 
to analyze the eyesight issues among the primary students considering 
several variables including gender, age, distance with devices, duration of 
using the devices, and analyzing previous eyesight issues. The findings 
indicated that male students have a more chance of developing eyesight 
issues than female students. Hence, it was deduced that female students who 
used their devices from a distance of more than 16 inches for less than 3 
hours had the least possibility of developing any eyesight issues. 


Keywords: COVID-19, eyesight issues, ophthalmology, logistic regression 
1. INTRODUCTION 


Throughout history, notable pandemics such as AIDS, dengue, cholera, 
plague, influenza, smallpox, and severe acute respiratory syndrome have 
had far-reaching consequences, affecting not only the healthcare system but 
also the economy, society, and security. Several key factors contribute to 
the facilitation of a pandemic, including its global reach, novelty, severity, 


* Corresponding Author: Basitali.bukc @ bahria.edu.pk 


School of Seience I — I UMT 55 


Volume 7 Issue 2, 2023 P 


Assessing the Impact of Prolonged... 


reduced immunity, high infection rates, and other related aspects [1]. 


Prior to early December 2019, numerous symptoms resembling 
pneumonia were documented in Wuhan, located in Hubei province, China 
[2]. These symptoms occurred because of the re-emergence of the 
coronavirus in the form of SARS-COV 2 [3]. In Pakistan, authorities 
reported the first two COVID-19 cases in March 2020 in Islamabad [4]. The 
lockdown was imposed countrywide and people were forced to isolate 
themselves in their houses. The traditional mode of education was shifted 
to an online education medium, which opened new gateways and new 
opportunities for students to learn this new medium of education [5]. Online 
education brought many challenges, which include a lack of proper 
infrastructure for educational setup, technical issues like internet lagging 
leading to voice and video disruption, and a lack of face-to-face interaction 
among teachers. Primary-level of students were dependent on their parents 
for their academic work and the online education system was new to their 
parents. Thereby, students or teachers developed several health conditions 
that were also connected to the online education system including eyestrain, 
cervical pain, lumbago, vertigo, headaches, psychological issues, lack of 
physical activity, and eyesight issues. 


The eye is the most delicate and sensitive organ of the human body. It 
is susceptible to any damage that could be chemical, physical, 
environmental, and others, which may affect its ability to envisage leading 
numerous eye-related issues. The structure of the eye mainly comprises the 
conjunctiva, sclera, iris, pupil, lens, and retina. The conjunctiva is a mucous 
membrane made up of non-keratinizing squamous epithelium containing 
goblet cells, langerhans cells, and dendritic melanocytes [6]. It protects and 
lubricates the eye by producing mucous and tears [7]. As light strikes the 
retina, several chemical reactions take place within different cell layers such 
as rods and cones. These chemical reactions in the form of electrical signals 
disseminate through nerve cells into the second cranial nerve [8]. Several 
abnormalities associated with the eye include ageing, increased exposure to 
unshielded screens, trauma, infections, and family history. The main issues 
that could occur due to prolonged exposure to screens are eye itching, 
irritation, burning, redness, drying, and fatigue. Long-sightedness 
(hypermetropia) and shortsightedness (myopia), which are the most 
common abnormalities in growing ages. Hypermetropia occurs when the 
eyeball is too short so the focus point is behind the retina at the back of the 
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eye. If this condition occurs in students, a squint may be noticed, for 
instance, one eye turning inward. On the other hand, myopia occurs when a 
person is unable to see only nearby objects.This condition occurs when the 
length of the eye is excessive, causing the light to converge in front of the 
retina located at the back of the eye. Myopia is a serious condition because 
if a child becomes myopic before nine years old, there are chances that he 
or she may develop a high level of myopia in the future, which makes them 
prone to glaucoma, retinal detachments, and myopic retinal degeneration. 
There are also chances of complete loss of eyesight for the high-level 
myopic person [9]. Dry eye is an eye condition that leads to ocular 
tenderness. It is characterized by a disruption in the tear film, resulting in 
the drying of tears. This condition is commonly associated with decreased 
intervals between blinking [10, 11]. Continuous use of screens can cause 
eye fatigue and soreness, which contributes to the eye injury and diseases 
like refractive error, myopia, cataract, and glaucoma [12]. Recently, eye 
abnormalities are associated with increased exposure time to screens among 
all age groups. In the last couple of years, as the COVID-19 lockdown 
prompted a global shift towards online education, numerous reports in the 
literature have highlighted the prevalence of eyesight problems. During the 
COVID-19 pandemic lockdown in the past two years, the traditional 
education system shifted to online educational medium. Reportedly, many 
eyesight cases were noted during this particular time of the COVID-19 
pandemic. 


However, in Pakistan, such issues remained unnoticed due to the dearth 
concerns, such as financial issues and regular checkup. People in this crucial 
time avoided visiting hospitals because of fear of getting COVID from the 
hospitals of infected patients. According to WHO, 2.2 billion people around 
the world suffered from eye conditions [13]. Bahkir reported that the 
increase in the usage of digital devices during COVID-19 led to an adverse 
effect on ocular health in all age groups. Therefore, 90.42% of the 
respondents reportedly experienced eye abnormalities due to the excessive 
usage of digital devices. In their study, they identified a high incidence of 
dry eye and headache in females. A percentage of 43.5% of people 
experienced headaches, 29% had eye pain, 23.8% had heavy eyelids, and 
74.2% had sleep disturbances [14]. Balsam Alabdulkader through a self- 
reported questionnaire estimated that the prolonged use of digital devices 
increased during the COVID-19 pandemic and reported that 78% of 
students faced digital eyestrains [15]. Bengi Demirayak elaborated increase 
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in digital device use during COVID-19 exacerbated the digital eyestrain 
problem in students. Hence, 52.2% of students experienced headaches as 
the most common symptom [16]. Abdullah Altalhi conducted an electronic 
self-administered survey on KSAUHS students in Saudi Arabia of 
Computer Vision Syndrome (CVS). It was reported that CVS symptoms 
were reportedly higher in female students who used to wear spectacles and 
observed glare on the screen. Reportedly, the most common symptoms 
noticed among students were headache (68%), feeling of affected eyesight 
(65%), eye itchiness (63%), burning sensation (62%), excessive blinking 
(40%), and difficulty moving eyelids (9%) [17]. 


In China, more than 180 million of primary and secondary school 
students were circumscribed to their homes. It was observed that students 
were less physically active with longer screen time and irregular diet 
resulted in cardiorespiratory weaknesses [18]. Homeschooling is a unique 
and alternative format of education. According to the student demographic 
questionnaire, 17.6% of students were identified as having behavioural and 
emotional effects. However, low-grade students were more at risk with 
19.3% as compared to high-grade students having 13.7% effects [19]. 
According to the Indian Journal of Ophthalmology, a survey conducted 
among 217 individuals showed the frequency of digital eyestrain to be 
50.23%. Common manifestations observed were itching in the eyes and 
headache (53.9%). This survey-based research highlighted the prevalence 
of digital eyestrain causing ocular health issues due to increased usage of 
digital devices [20]. Increased usage of digital devices among students for 
online sessions showed the progression of myopia in individuals. The 
percentage of myopic individuals reached up to 10.40% in Chongqing, 
China [21]. Researchers used the chi-square test for analyzing the immense 
use of computers by university students and it was identified that burning 
eyes and tired eyes were the common symptoms reported among students 
[22]. Computer vision syndrome was reported in 75% of visual display 
terminal workers as compared to other workers [23]. 


Ugam Usgaonkar through his e-survey-based research estimated that 
most of the respondents were involved in more than two digital activities on 
more than two digital devices. The symptoms observed were watery eyes, 
headache, dry eyes, shoulder pain, pain behind the eye, and back pain. The 
increase in ocular symptoms was due to an increase in time spent on digital 
devices [24]. Richa Gupta elaborated on the adverse effect of e-learning on 
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the eye health of school students during COVID-19 through survey-based 
research. A percentage of 92.8% of students reportedly had eyestrain 
symptoms. The most common symptoms were eye redness, which was 
69.1% and heaviness of the eyelid, which was 79.7% [25]. Pratyusha Ganne 
analyzed the digital eyestrain problem due to the increase in digital device 
usage during COVID-19. The prevalence of eyestrain was found to be 
higher in students, which were attending online classes. The prime factors, 
which were involved, were pre-existing eye disease, greater screen time, 
screen distance, and using the screen in the dark [26]. 


In this paper, a survey-based research regarding the ophthalmic health 
of primary school students was conducted as a result of attending online 
classes through different digital devices during COVID-19. The current 
study was designed to analyze the impact of online education on their 
eyesight is qualitative, based on survey of 100 primary students, 59 female, 
and 41 male students conducted for the age group 4 to 12 years. It was also 
observed that some students were already facing eye abnormalities, which 
may include eyestrain, red eyes, night blindness, uveitis, presbyopia, 
floaters, dry eyes, cataracts, retinal disorders, corneal diseases, and others. 
For this purpose, Logistic regression analysis would help to estimate the 
probability of the eyesight issues of the primary students after online 
education. This dichotomous variable would be used as a dependent 
variable, which depends on gender, age, the distance of the device with the 
eyes, duration of the user device, and already facing any eye health issues. 


2. METHODOLOGY 
2.1. Statistical Model and Analysis 


The relationship between a dichotomous dependent variable and one or 
more independent categorical variables was illustrated by binary logistic 
regression analysis. This model provides odds, odds ratio, and probability 
Statistics, which can be highly useful in properly expressing the research 
goal [27]. Generally, this model is used while discussing a dependent 
variable, which is binary, dichotomous, or alternative, rather than 
continuous variables [28, 29]. There are other statistical models, which can 
be applied in previous literature related to statistical convergence [30, 31]. 


The IBM SPSS version 25.0 software was used to find out the 
relationship between eyesight issues with the explanatory variables of 
gender, age, the distance of the device with the eyes, duration of using the 
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device, and people who were already facing any eye abnormalities. In the 
study population 100 primary students who have taken online classes during 
the pandemic, 59 were female and 41 were male students with different age 
groups ranging from 4-12 years. The device distance could be less or greater 
than 16 inches. The duration of using devices for online education was 
noticeably more or less than 3 hours. Another significant aspect to consider 
is whether the student had experienced any eye-related issues prior to the 
onset of the COVID-19 pandemic. After executing the logistic regression in 
IBM SPSS, some information was obtained, which helped to estimate the 
impact of the devices used for online education on primary students’ eyes. 


2.2. Model Diagnosis 


Model diagnosis shows the complete analysis of the model execution. 
Table 1 provides a complete indication of how well the model executes. 
This is also termed as а "Goodness of fit test". It 
is utilized to decide the rate of an exact number 
of expectations in total predictions. The significant value in the Omnibus 
test is less than 0.05, which reflects that the model is significant. The chi- 
square statistic is 19.224 with 5 degrees of freedom. Therefore, the overall 
model is statistically significant, 72 (5) = 19.224, p < 0.05. The Hosmer and 
Lemeshow tests were also used to support the current model but their 
interpretation is completely different from the omnibus test. In this test, the 
significance value should be higher than 0.05 so that the model does fit the 
data. Table 1 shows that the chi-square value for the Hosmer and Lemeshow 
Test is 11.980 with a significance level of 0.101. This value is larger than 
0.05; therefore, indicating accuracy for the model. 


Table 1. Omnibus & Hosmer Tests of Model Coefficients 


Test Chi-square Degrees of Significant 
freedom value 
Omnibus 19.224 5 0.002 
Hosmer and 11.980 2 Me 
Lemeshow 


To calculate the variations in the dependent variable, the method of the 
Cox & Snell R? and Nagelkerke R? were supportive. According to these 
methods, Table 2 shows that around 17.5%-23.5% of the variability in the 
eyesight issues after the COVID-19 pandemic were considered by all the 
independent variables. 
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Table 2. Model Summary 


-2 Log likelihood Cox & Snell R Square Nagelkerke R Square 
116.835 0.175 0.235 


2.3. Model Characteristics 


The classification table depicted in Table 3, which is obtained after the 
execution of logistic regression in IBM SPSS, explains how well the model 
is going to predict the correct category (facing eyesight issues/no eyesight 
issues after the COVID-19 pandemic) for each case. The sensitivity of the 
model is the percentage of the students having eyesight issues after COVID- 
19 pandemic, which were identified by the model in Table 3, which is 
primarily 52.496 of the students who are facing eyesight issues after 
COVID-19 pandemic. 


Table 3. Classification Table 


Predicted 
Is the child facing any 
Observed eyesight issues after Percentage 
COVID-19 pandemic? Correct 
no yes 
Is the child facing any no 44 14 75.9 
eyesight issues after 
COVID-19 pandemic? e a x Би 
Overall Percentage 66.0 


The model's specificity refers to the percentage of students who do not 
have eyesight issues after COVID-19 pandemic. In this model, the 
specificity is 75.9%. The model accurately predicted that students who did 
not have eyesight problems after the COVID-19 pandemic were indeed free 
of such issues. 


3. RESULTS AND DISCUSSION 


The model shows a complete picture of all the parameters measured in 
this study. Table 4 provides a detaileddata about the contributions of every 
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independent variable. 


Table 4. Coefficients of Logistic Regression Model for Different Eyesight 
Issues Parameters 


95% C.I. for 
Independent Coefficients 
Б SE. Wald df Sig. Ехр(В) EXP(B) 
variables (B) 
Lower Upper 
Gender of child 0.774 0.465 2.774 1 0.096 2.168 872 5.389 
Age of child 0.204 0.107 3.616 1 0.057 1.226 ‚994 1.513 


What was the 

distance of eyes of 

the child from -0.748 0.581 1.653 1 0199 0474 0.152 1.480 
computer/phone 


screen? 


What is the 
maximum hours of 

0.876 0.517 2.867 1 0.090 2.402 0.871 6.622 
screen exposure of 


the child? 


Did the child have 
any eyesight issue 


before COVID-19 


0.903 0.484 3.478 1 0.062 2466 0.955 6.369 


pandemic? 


Constant -3.206 0.998 10.316 1 0.001 0.041 


The direction and quality of the relationship between the dependent and 
independent variables are shown in the column B. In the model, the B value 
of the distance of eyes of the students from computer/phone screen is 
negative, which means that the distance is more than 16 inches and has less 
probability of having eyesight issues than if the distance is less than 16 
inches. The other B values are positive, which reflects that if the age of 
students increased, the chances of having eyesight issues increases. 
Furthermore, the B value of gender is also positive, which means that 
female students has less probability of having eyesight issues than male 
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students. More than 3 hours of screening time is more harmful for the 
students and create more eyesight issues than the other students whose 
screening time is less than 3 hours. Similarly, students who were already 
facing vision problems have more possibility to have more severe vision 
problems after online education. 


These values of B can be utilized in an equation to find the probability 
of the primary students having eyesight issues. These values are in log- 
odds units, which can be converted to probabilities. The following is the 
expression for the binary logistic regression model considering the 
aforementioned variables: 


log С) = bo + by Xgenaer + 02 Харе + D3 Xaistance + 04 Xduration + 

bs Xalready eyesight issues (3.1) 
Where p, is the probability of having eyesight issues after COVID-19 
pandemic log G5) in Eq. (3.1) is considered as a log of odds or link logit 
function while rm is taken as the odds ratio. The logistic regression 


equation in variables form is: 


log (+) —- 3.206 + 0.774 Xpender + 0.204 Xage — 0. 748Xaistance 


+ 0.876 Xduration + 0.903 Xalready eye issues (3.2) 
Applying antilog on both sides, the odds ratio will be equal to, 
[= 
1-р 
= е-3-206+0.774 Xgender + 0.204 харе —0.748X distance #0.876 Xquration + 0.903 Xalready eye issues 
(3.3) 
The probability of each eyesight issue parameter with different 


covariates can be easily depicted. Eq. (3.4) represents the probability 
equation generated from the model. 


e 3.206+0.774 Xgender + 0.204 Xage—0.748Xdistance + 0.876 Xauration * 0.903 Xalready eyesight issues 
E 1+ e 3:206--0.774 Xgender * 0.204 Xage—0.748Xdistance + 0.876 Xquration f 9-903 Xalready eyesight issues ` ` 


(3.4) 
In Table 4, column Exp (B) is the value of the odd ratio for each 
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independent variable. It shows the increase in probability of occurrence as 
the independent variable increases by one unit. In the model, the odds of a 
student having eyesight issues is 0.474 times higher for the student whose 
distance of eyes from the device is less than 16 inches than a child whose 
distance of eyes from the device is more than 16 inches, if all other factors 
are supposed to be constant. Similarly, the odds of a student having eyesight 
issues is 2.402 times higher for the student whose screen time is more than 
3 hours than those students having screen time is less than 3 hours. Same as 
in the situation of suffering any other eyesight issues before COVID-19 
pandemic, those students who are suffering any eyesight issues before 
COVID-19 pandemic, the eyesight issues is 2.466 times higher than those 
who not suffering any eyesight issues before online education, taking all 
other factors constant. The odd ratio of the gender is greater than 1, which 
represents that mostly male students have reported eyesight issues than 
female students. The odds of eyesight issues in female students decrease 
by a factor of 2.168 to male students when all other factors are kept constant. 
In the model, one continuous variable, age, is also included. The odd ratio 
of age is increased by a factor of 1.226 for a one-year increase in age. In 
Table 4, the 95% confidence intervals provide lower and upper bounds 
associated with the corresponding odds ratios. These are the limits where 
you can be 95% confident that they cover the true value of the odds ratio. 
Using the probability equation obtained in (4), Table 5 shows that all 
probabilities of having eyesight issues were measured with the impact of all 
independent variables. 


As age is the continuous variable, maximum, minimum, and average 
age has been taken to estimate the probabilities. The lowest probability of 
having eyesight issues is 0.04156, which depicts a female primary student 
who is not suffering any eyesight issues before online education, at the age 
of 4 years, with the distance of more than 16 inches with screen and her 
exposure duration is less than 3 hours. While the maximum probability of 
having eyesight issues is 0.85753, which represents a male student who is 
already facing eyesight issues at the age of approximately 12 years, with a 
distance of fewer than 16 inches and his exposure duration is more than 3 
hours. Figure 1 shows the probability of having eyesight issues in male 
students with the consideration of all the independent variables. While, 
Figure 2 shows the probability of having eyesight issues in female students 
with the consideration of all the independent variables. 
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Table 5. Probability Table of Having Eyesight Issues with the Covariates 


Distance Less than 16 inches More than 16 inches 


Duration Less than 3 hours More than 3 hours Less than 3 hours More than 3 hours 


Eyesight condition 


before pandemic Healthy Unhealthy Healthy Unhealthy Healthy Unhealthy Healthy Unhealthy 

4 year 0.083 0.184 0.180 0.351 0.041 0.096 0.094 0.204 

E <8 year 0.171 0.338 0.332 0.551 0.089 0.194 0.190 0.367 
m 

X12 year 0.319 0.536 0.529 0.735 0.181 0.353 0.347 0.567 

4 year 0.165 0.328 0.323 0.540 0.085 0.188 0.184 0.357 

3 <8 year 0.310 0.525 0.518 0.726 0.175 0.344 0.338 0.557 

<12 year 0.503 0.714 0.709 0.857 0.324 0.542 0.535 0.740 
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m4 year Ш<8 year m<12 year 


0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


Healthy Unhealthy Healthy Unhealthy Healthy Unhealthy Healthy Unhealthy 
eyes eyes eyes eyes eyes eyes eyes eyes 


Probability 


Less than 3 hours | Greater than 3 hours | Less than 3 hours Greater than 3 hours 


Figure 1. Probability of Eyesight Issues in Male Students 


m4year Ш<8 year 


Healthy Unhealthy] Healthy Unhealthy] Healthy Unhealthy} Healthy Unhealthy 
eyes eyes eyes eyes eyes eyes eyes eyes 


Less than 3 hours | Greater than 3 hours | Less than 3 hours | Greater than 3 hours 


Figure 2. Probability of Eyesight Issues in Female Students 


3.1. Conclusion 


The current study deployed a logistic regression analysis, which 
estimated the probability of optical impairment in primary school children. 
Noticeably, these children were taking online classes during the COVID-19 
pandemic. This research paper indicated that a logistic regression model 
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considered the explanatory variables of gender, age, eye distance from the 
device, duration of use, and previous vision problems to model eyesight 
abnormalities in primary school children. With the help of this model, it has 
been observed that male students have more chances than female students 
of developing eyesight issues. Furthermore, it was observed that the 
probability of experiencing eye problems increases when factors such as 
age, distance, and duration of online classes are extended. Additionally, if 
the child has preexisting eye disorders, the likelihood of developing eye 
issues increase. Thereby, it was identified that female students of 4 years of 
age have the lowest chance of having eyesight problems if they use devices 
with 16 inches of distance. This is because they are exposed for less than 3 
hours and do not have any prior eye problems. Male students ranging 
between 9-12 years, using a screen with a distance of less than 16 inches 
over an extended period of up to 3 hours and who have already faced vision 
issues have the highest risk of having vision issues. Therefore, it is 
recommended that male students should wear eyeglasses or should maintain 
a proper screen distance to reduce the impact on their eyesight. 
Additionally, periodic eye exams enable them to detect earlier symptoms of 
any developing problems. 
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ABSTRACT 


Lawsone is an orange red dye mainly isolated from Lawsonia inermis plant 
leaves. L. inermis is commonly known as henna plant, Egyptian privet, or 
mignonette tree. In Pakistan, henna is known as mehndi. Paste from thick 
henna leaves has been used to dye hairs, fingernails, palms of hands, and 
eyebrows for more than 5000 years. The characteristic colour of henna is 
because of the natural dye lawsone present in it. In this study, lawsone was 
extracted from four samples of henna leaves collected from various parts of 
Pakistan. Moreover, the impact of the environment on the leaves was 
examined. Lawsone was extracted using diethyl ether, so that the maximum 
amount of lawsone could be isolated from the leaves. The maximum and 
minimum amount of lawsone isolated during this work was 1.01% and 
0.745% of dry weight, respectively. Afterwards, IR and HPLC were used to 
identify the isolated reddish-brown substance known as lawsone which is 
useful in textile, cosmetics, and pharmaceutical sectors. 


Keywords: dye, Henna, Lawsonia inermis, lawsone, plant 
1. INTRODUCTION 


Lawsonia inermis is commonly known as henna plant, Egyptian privet, 
or mignonette tree [1, 2]. In Pakistan, henna is known as mehndi. The 
principal colouring component of this plant is known as lawsone. This plant 
is grown on a large scale in Sindh and Baluchistan, two provinces of 
Pakistan [3]. However, due to its medicinal and cosmetic applications, it is 
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cultivated throughout the world [4]. It is a significant medicinal plant that 
has drawn the attention of scientists from all over the world due to its 
amazing pharmacological properties [5]. It is a mid-sized, highly branched 
shrub. Several medicinal uses have been reported for its flower, seed, roots, 
and stem [6]. It is used to dye hair, skin, and nails [7]. Skin and hair colour 
are extracted from lawsone by using the technique based on Michael 
addition reaction theory. A non-permanent colouration that can last up to 
two to three weeks is produced when Lawsone makes a covalent link with 
the keratin protein present in skin and hair. The interaction between lawsone 
and keratin determines colour intensity [8]. Darker colouration results from 
a greater lawsone-keratin interaction. However, with the passage of time, 
lawsone is dissolved in water. As a result, lawsone concentration decreases 
and henna stain fades away [9]. 


Chemically, lawsone is a 2-hydroxy-1,4-naphthoquinone, also known 
as hennotannic acid [10]. Lawsone is an orange-red dye mainly found in 
Lawsonia inermis leaves. It is also present in other plant components, 
although the amount is quite small [11]. Lawsone gets its name from the 
plant it is derived from, that is, Lawsonia inermis. The genus Lawsonia 
includes only two species, that is, L. inermis and L. odorata. Lawsone is 
present in very small amounts in Lawsonia odorata. A minor quantity of 
lawsone was also observed in water hyacinth bloom [12]. 


Human beings have been using henna extracts containing lawsone for 
at least 5000 years to dye their hair and skin [13]. Lawsone shows an 
excellent antifungal, antibacterial, antitumor, antiviral, and pesticidal 
activity [14, 15]. The amount of lawsone in henna plant leaves varies by 
region. Geological conditions and a number of environmental factors affect 
its concentration [16]. Moreover, soil temperatures between 25°C and 30°C 
are ideal for the germination of this plant. Although, it may also survive in 
the temperature range of 35°C to 45°C [17]. In comparison to winter, a 
higher concentration of lawsone is observed in henna leaves during summer 
[18]. The current study highlights the effect of environmental conditions on 
lawsone content cultivated in different areas of Pakistan. 


2. MATERIALS AND METHOD 


Experimental work regarding the isolation and identification of lawsone 
content from the leaves of henna was carried out in the laboratory of 
Pakistan Council of Scientific and Industrial Research (PCSIR), Lahore. 
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HPLC instrument used in this experiment was designed by Shimadzu 
Corporation, Japan. While, the column used was LC-18 and SPD 6AV UV- 
detector, adjusted at 260nm. Henna leaves used in this work were collected 
from different areas of Sindh and Punjab. The identification of samples was 
confirmed by the Department of Botany, University of Education, Lahore. 
All the samples were shade dried. 


Dried henna leaves sample (100g) was added in a 5L Erlenmeyer flask 
containing 5000ml distilled water and refluxed at 80°C till the vapours 
stopped coming out of the flask. The suspension was stirred continuously 
with the help of magnetic stirrer to avoid any sort of coagulation. After 
approximately 3 hours, the colour of suspension turned brown and the reflux 
assembly was removed. Suspension was left overnight. Then, the 
suspension was made alkaline by adding sodium bicarbonate (20g). 
Alkaline suspension was then filtered and the filtrate was made acidic by 
adding 0.12M hydrochloric acid till the pH of the solution became 3. 
Lawsone was extracted by using diethyl ether. Excessive solvent was 
evaporated on a hot plate at 30°C to avoid the burning of lawsone. Isolated 
reddish-brown material was then scratched and stored in a pre-weigh vial. 
The same procedure was reported for other samples and the isolated 
lawsone was stored in pre-weighed vials, separately (Figure 1 and 2) [16]. 
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3. RESULTS AND DISCUSSION 


The amount of lawsone isolated from four different henna leaves 
samples was not more than 1% of the dry weight of samples (Table 1). The 
percentage of the isolated reddish-brown material depicted the effect of 
environmental conditions on the lawsone content. The material was 
identified by using IR and HPLC and the results were compared with the 
standard chemical lawsone. 


Table 1. Percentage of Lawsone Isolated from Different Henna Samples 


Henna Leaves Samples Percentage of Lawsone Extracted 
Sample leaves | (S1) 1.010 
Sample leaves 2 (S2) 0.845 
Sample leaves 3 (S3) 0.785 
Sample leaves 4 (S4) 0.744 


Figure 2. Reddish Brown Material Isolated from Different Henna Leaves 
Samples 
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3.1. IR Spectra 


The IR spectra of the isolated material shows different peaks. Each peak 
corresponds to a specific functional group. The stretching of the O-H 
functional group, overlaid by the naphthalene ring’s C-H vibrations, 
reached its maximum at 3159 cm'!. Carbonyl stretching caused реак 
formation at 1675 cm-1 and 1636 cm'!. These peaks were expected to divide 
due to internal hydrogen bonds. Peaks at 1216 cm’! were related to the 
stretching C-O, while peaks at 1584 and 1510 сп! were related to the 
naphthalene ring's C=C vibrational bands. The IR spectra peaks confirmed 
the chemical nature of the reddish-brown substance, namely lawsone 
(Figure 3). 
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Figure 3. IR Spectra of Isolated Reddish-Brown Material 
3.2. HPLC Analysis 


The nature of the isolated material was further confirmed by HPLC. For 
this purpose, standard lawsone was purchased from Merck. A total of 50mg 
of standard and isolated material was used to prepare standard and sample 
solutions, respectively. Mobile phase used during the experiment was 70% 
methanol and 30% water. A significant resemblance between the separated 
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material’s HPLC spectra and the standard lawsone HPLC spectra was 
observed. This resemblance proved that the isolated chemical was lawsone 
(Figure 4). 


Figure 4. HPLC Spectra of Standard and Isolated Reddish-Brown 
Material 


3.3. Effect of Geological Conditions on Lawsone Content 


Henna leaves collected from areas with a high temperature and dry air 
have more lawsone. High temperature is favorable for lawsone content in 
henna leaves. Moreover, henna plants cultivated during summer have more 
lawsone content in their leaves than those cultivated in winter. 


3.4. Factors Effecting Lawsone Content 


3.4.1. Effect of Temperature. Temperature has a strong influence on 
the amount of lawsone extracted from Lawsonia inermis leaves [19]. This 
is because if the extraction temperature exceeds 80°C, then most of lawsone 
is either dissolved in water or burned out. Similarly, if the temperature 
exceeds the boiling temperature of the solvent during the isolation of 
lawsone from diethyl ether, then all the lawsone is burnt. 


3.4.2. Effect of pH. Change in pH does affect lawsone concentration 
but its effects are less harsh as compared to temperature. A slight change in 
pH can be neglected but a major change does affect the lawsone yield. The 
greater the pH of solution the lesser is the amount of the isolated lawsone. 
pH also affects the colour of the extract [20]. 
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4. CONCLUSION 


The primary colouring component of the Lawsonia inermis plant is 
lawsone. The identification of environmental conditions and their impact on 
lawsone content were studied in this research. Lawsone was extracted from 
four different kinds of henna leaves samples. The maximum and minimum 
amount of lawsone isolated during this work was 1.01% and 0.745% of dry 
weight, respectively. The percentage yield made it evident that sample 
leaves from Sindh contained more lawsone than sample leaves from Punjab. 
Thus, it was concluded that the weather conditions of Sindh are in favor of 
lawsone content in Lawsonia inermis leaves. Moreover, it was also 
determined that henna plants harvested in the summer season have more 
lawsone than those harvested in winters. Thus, isolated lawsone can be used 
in the pharmaceutical industry and as a dying agent in cosmetic and textile 
industries. 
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ABSTRACT 


Ephedra gerardiana (E. gerardiana,), commonly known as "Ma-huang," is 
a medicinal plant native to the Himalayan region. Historically, it has a long 
traditional use in various indigenous systems for medicinal purpose. The 
current review aims to summarize the present knowledge on the botanical 
description, phytochemistry, traditional uses, and pharmacological 
properties of E. gerardiana. 'This study was thoroughly based on previous 
literature and was conducted by using electronic databases such as PubMed, 
Scopus, and Google Scholar. Relevant studies and recently published 
articles up till 2020 were included in this review. E. gerardiana is a small 
shrub-like plant characterized by jointed stems, scale-like leaves, and 
cylindrical cones. It belongs to the family Ephedraceae and is commonly 
found in arid and semi-arid regions. Phytochemical analysis has revealed 
the presence of various bioactive compounds in Ephedra gerardiana, 
including ephedrine, pseudoephedrine, norephedrine, and methyl 
ephedrine. These alkaloids are primarily responsible for their 
pharmacological activities. In traditional medicine, Ephedra gerardiana has 
been used for its bronchodilatory, diaphoretic, diuretic, and thermogenic 
properties. Additionally, it has been employed for the treatment of 
respiratory disorders, namely, asthma and nasal congestion, as a stimulant. 
Its bioactive alkaloids, particularly ephedrine, have been extensively 
studied for their potential therapeutic applications. However, caution should 
be exercised due to the potential side effects associated with ephedrine- 
containing products. Further research can explore its mechanisms of action, 
optimize dosage regimens, and investigate potential drug interactions to 
ensure its safe and effective use in clinical settings. 
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1. INTRODUCTION 


Ephedra gerardiana is a herbal drug, which contains many useful 
features. This herbal drug has many medicinal uses. It belongs to the family 
of Ephedraceae and its plant is approximately 1-4 feet long. It has red 
berries because of this; the plant becomes more beautiful and attractive. 


The red colour of berries is due to the presence of their chemical 
constituent Flavonones. 


According to the Unani system of medicine temperament of ephedra is 
hot and dry in a second order. In the Unani system of medicine drugs, which 
belong to second-order temperament, are used to treat phlegmatic disorder 
more easily and rapidly. 


The pharmacological activities include anti-inflammatory, alterative, 
diuretic, expectorant, anti-inflammatory, and bronchitis. 


The most effective action of this herbal drug is on our respiratory 
system. In Chinese medicine, it is also known as Mahaung. 


Its most important biomarkers are Ephedrine and Pseudo-ephedrine. 


The plant parts, which are used in medicine, are branches, flowers, and 
fruits. 


2. METHODOLOGY 


Different online databases were examined to highlight the therapeutical 
potential of Ephedra gerardiana, including Springer-Link, Google, 
PubMed, MEDLINE Data, Wiley Online Library, Science Direct, and 
Google Scholar. Data were searched using different keywords, such as 
pharmacological/biological activities, chemical constituents, morphology, 
habitat, characteristics, and “distribution of Ephedra gerardiana. 


2.1. Pharmacognostic Features 


Latin Name: Epheragerardiana Wall, Ephedra vulgaris 
Family: Epheraceae. 

Arabic Name: ZanbulKhail. 

Urdu Name: Somkalpana, Masania, Narom, Somkalpa 
English Name: Ephedra, Joint Pine, Brigham Tea. 


School of Sieme | —3à$i UMT— 83 


Volume 7 Issue 2, 2023 ae 


Review of Pharmacological Activity... 


Hindi Name: Som, Tutgantha. 


2.2. Characteristics 


It is a small shrub, which grows up to 1-4 feet. The stem of the plant has 
a nodule with dark green branches. 


The leaves are small and round. The flowers are 6-10 mm long and 
unisexual naturally. The fruit of the plant is sweet, round in shape, red in 
colour, and has 1-2 seeds. 


2.3. Habitat 


Endemic to the mountains of Afghanistan, Bhutan, Northem India, 
Nepal, Sikkim, and Tibet. 


2.4. Chemical Constituent 


In all types of herbal drugs, there are different chemicals, which give 
them different pharmacological activity and medicinal uses, along with 
different side effects. 


The ephedra chemicals are: 
lignans, flavonoids, pseudoephedrine, ephedrine. 


The root contains spermine alkaloids, Phenylethylamine alkaloids, 
Polysaccharides, Catechols, Tyrosine Derivatives, and EphedranineA& B. 


The most important chemicals are ephedrine and pseudo-ephedrine [1— 
3]. 


The chemical constituents , which are separated from the herbal plants 
are more than 145 including alkaloids, flavonoids, tannins, and 
polysaccharides. The compounds separated from each species are enlisted 
in Tables 1-4. Their chemical structures are given in Fig. 1—5. 


Table 1. Alkaloids Different Species of the Genus Ephedra (the structures 
of main compounds are illustrated in Fig. 1) 


Classification Chemical name Plant source Ref 
Macrocyclic Ephedradine A (1) E. sinica (root) [4] 
spermine Ephedradine B (2) E. sinica (root) [5] 
alkaloids Ephedradine C (3) E. sinica (root) [5] 
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Classification Chemical name Plant source Ref 
Ephedradine D (4) E. sinica (root) 
Imidazole | : w: 
alkaloid Feruloylhistamine (5) E. sinica (root) [6] 
D(-)-Ephedrine (6) 
L(+)-Pseudoephedrine (7) E.Sinica,E.nebrodensis, 
D(-)Norephedrine (8) Fe major [7] 
: E. sinica, E. nebrodensis ~~ 
L(+)-Noreseudoephedrine (9) un i . [8] 
D(-)Methyl ephedrine (10) E. udi E. nebrodensis [8] 
L(+)-Methylpseudoephedrine(11) и [9] 
Ephedroxane (12) Е. sinica, E. nebrodensis < 
Атр hetamine- : m E. sinica, E. nebrodensis [8] 
type alkaloids 3, 4-Dimethyl-5-pheyloxazolidine [10] 
(13) E. sinica "T 
2, 3, 4-Trimethyl-5- E. sinica m 
phenyloxazolidine (14) E. sinica A 
O-benzoyl-L(+)-pseudoephedrine E. sinica ІШІ 
(15) E. sinica ШІ 
O-benzoyl-D(—)-ephedrine (16) E. aphylla 
Hordenine (17) 
E. transitoria 
Transtorine (18) E. pachycladassp. [12] 
6-Methoxykynurenic acid (19) Sinaica [13] 
Quinoline : i 
: Kynurenic acid (20) E. pachycladassp. [13] 
alkaloids DP 
6-Hydroxykynurenic acids (21) Sinaica [13] 
Ephedralone (22) E. foemineassp. Foliate [14] 
E. Alata, E. aphylla 
Pyrrolidine ; : E. foemineassp. foliate 
alkaloid cis-3, 4-Methanoproline (23) (seed) [13] 
Maokonine (24) 
(+)-1-Phenyl-2-imido-I-propanol Ё. sinica(root) [15] 
(25) E. sinica [16] 
Other ET 5-4 
alkaloids Tetramethylpyrazine (26) E. sinica [16] 
N-methybenzlamine (27) E. sinica [16] 
(2,5,35,45)- Е. altissima [13] 


2(carboxycyclopropyl)glycine (28) 
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Figure 1. Structures of Different Alkaloids of Ephedra 
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Table 2. Flavonoids of Different Species of the Genus Ephedra (the 
structures of main compounds are illustrated in Fig. 2) 


Classification Chemical name Plant source Ref. 
Herbacetin (29) 
Kaempferol (30) E. sinica 
Quercetin (31) E. 
Herbacetin 7-methyl ether (32) sinica(root) [17] 
Rutin (33) E. [18] 
Herbacetin 8-methyl ether3-O- sinica(root) — [18] 
glucoside-7-O-rutinoside (34) E. sinica [17] 
Herbacetin 7-O-(6"-quinylglucoside) Е. [9] 
(35) campylopoda [19] 
Herbacetin 3-O-rhamnoside 8-O- E. alata [19] 
glucoside (36) E alata [19] 
Flavonols Pollenitin B (37) E. aphylla [20] 
Herbacetin-8-methyl ether 3-O- E. sinica 
lucosid | [20] 
glucoside (38) ЖЕ 
i : E. sinica [20] 
Herbacetin 7-O-glucoside (39) - 
. E. sinica [20] 
Kaempferol 3-O-rhamnoside 7-O- Е sini 
glucoside (40) то [20] 
Herbacetin 7-O-neohesperidoside Е сак [20] 
(41) Е. sinica [19] 
Kaempferol-3-O-glucoside-7-O- E. alata [19] 
rhamnoside (42) E. alata [9] 
Kaempferol 3-O-rhamnoside (43) E. 
Quercetin 3-O-rhamnoside (44) campylopoda 
Quercetin-3-O-glucoside (45) 
E. 
. Dihydroquercetin (46) sinica(root) [18] 
Dihydroflavonol . { 
3-Hydroxynaringenin (47) E. [18] 
sinica(root) 
3’, 4’, 5, 7-Tetrahydroxy flavanone E. 
° ° ЫЛ j 8 
гі (48) sinica(root) = 
е Naringenin (49) E. [21] 
Hesperidin (50) sinica(root) 1221 
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Classification Chemical name Plant source Ref. 
E. sinica 
Flavanols (—)-epicatechin (51) E. sinica E. [23] 
(—-epiafzelechin (52) nebrodensis [18] 
Catechin (53) E. | [20] 
Afzelechin (54) pos dn [18] 
Leucocyanidin (55) В а [20] 
Symplocoside (56) sinicatroub [20] 
E. sinica 
E. sinica 
г 5; Е. sinica 
Tricin (57) 
Luteolin (58) alee san 
Apigeni (59 F ыш 
pigeni (59) | sinica(root) [18] 
COND Cu ae (60) БЕ Sinica [24] 
nb nio nas (61) E sinica [24] 
" Wu S PUDE SOHO) (62) E sinica [26] 
avones 
Жын. | Е. sinica [27] 
Неше Слао (64) E sien [20] 
тантана ы (65) Е [9] 
Vitexin (66) campylopoda [50] 
Lucenin III (67) E. sinica [19] 
2", 2"-Di-O-f-glucopyranosyl- E. alata [28] 
vicenin II (68) 
E. aphylla 
mum Leucodelphinidin (69) E. sinica [29] 
u Leucopelargonin(70) E. sinica [29] 
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OH 52 R,- R,- OH 54 R, - R,- OH 
53 R, =H R,- OH 55 R, =H R,=OH 


OH 
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56 59 R,=OH R,-H a 
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Figure 2. Structure of Different Flavonoids of Ephedra 
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Table 3. Tannins of Different Species of the Genus Ephedra 


Classification ChemicalName PlantSource Ref 
Ephedrannin A(71) 
Ephedrannin B (72) 
Мишай a) вабо: (al 
Muhuannin D (74) и 
à E. sinica(root) [30] 
Muhuannin B (75) 1 
à E. sinica(root) [31] 
Muhuannin E (76) da 
. E. sinica(root) [30] 
Muhuannin C (77) 2s 
: E. sinica(root) [30] 
Muhuannin F (78) n 
i E. sinica(root) [30] 
Muhuannin G (79) 5; 
: E. sinica(root) [30] 
Muhuannin H (80) чи 
š E. sinica(root) [32] 
Muhuannin I (81) Е 
: E. sinica(root) [32] 
Muhuannin J (82) a 
. E. sinica(root) [22] 
Muhuannin K (83) 
à E. sinica(root) [32] 
Ephedrannin D1 (84) ae 
. E. sinica(root) [32] 
Ephedrannin D2 (85) AM 
. А Е. sinica(root) [32] 
Dimer ргоап- Ephedrannin D3 (86) ч 
re . E. sinica [33] 
thocyanidins Ephedrannin D4 (87) i 
f E. sinica [33] 
Ephedrannin D5 (88) ie 
: E. sinica [33] 
Ephedrannin D6 (89) du 
2 Е. віпіса [33] 
Ephedrannin D7 (90) а 
š E. sinica [33] 
Ephedrannin D8 (91) 2 
2 Е. віпіса [33] 
Ephedrannin D9 (92) s 
: E. sinica [33] 
Ephedrannin D10 ЖЕ 
Е. віпіса [33] 
(93) А Е. sinica [33] 
Ephedrannin D11 i p 
E. sinica [33] 
(94) : E. sinica [33] 
Ephedrannin D12 СА a^ 
E. sinica [33] 
equ D13 E Sieg БЕЛ 
(96) Е. herb [33] 
Ephedrannin D14 
(97) 
Ephedrannin Tr1 (98) Е. sinica [33] 
Trimer proantho- Ephedrannin Tr2 (99) Е. sinica [33] 
cyanidins Ephedrannin Tr3 E. sinica [33] 
(100) E. sinica [33] 
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Classification ChemicalName PlantSource Ref 
Ephedrannin Tr4 E. sinica [33] 
(101) E. sinica [33] 
Ephedrannin Tr5 E. sinica [33] 
(102) E. sinica [33] 
Ephedrannin Tr6 E. sinica [33] 
(103) E. sinica [33] 
Ephedrannin Tr7 E. sinica [33] 
(104) E. sinica [33] 
Ephedrannin Tr8 E. sinica [33] 
(105) E. sinica [33] 
Ephedrannin Tr9 E. sinica [33] 
(106) 

Ephedrannin Tr10 

(107) 

Ephedrannin Тг11 

(108) 

Ephedrannin Tr12 

(109) 

Ephedrannin Tr13 

(110) 

Ephedrannin Tr14 

(111) 

Ephedrannin Tr15 

(112) 

Ephedrannin Tel 

(113) 

Ephedrannin Te2 ee [20] 

(114) uA 

А Е. sinica [20] 
Tetramer Ephedrannin Te3 i 20] 
proanthocya- nidins (115) S rose [ag 
: E. sinica [20] 

Ephedrannin Te4 E. sinica [20] 

(116) | == 

Ephedrannin Те5 

(117) 


Scientific Inquiry and Review 


2—SIR 


Volume 7 Issue 2, 2023 


76 R=OH 
77R=H 


92 R, = OH R, 
93 R, = OH R,= OH 


94R,-H 


ÓH 
о он 
HO. К 
ÓH 
HO 
96 R, =H 99 Rj-H R= 
97 R, = 0H 100 R,-OH R,=H 


School of Science 
Volume 7 Issue 2, 2023 


Review of Pharmacological Activity... 


104R,-H R,-H 
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Figure 3. Different Tannins from Ephedra 
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3. PHARMACOLOGICAL ACTION 
3.1. Cardiac Tonic 


Ephedra is used as a cardiac tonic. It affects the CVS due to its 
biomarker Ephedrine, which is sympathomimetic amine activity. It 
increase the heart rate and also helps to manage low blood pressure to treat 
bradycardia, concerning hypotension. It also increases cardiac contractility 
mostly in hypotension conditions in which blood pressure usually droops 


[1]. 


A case of orthotic hypotension was observed in this condition having 
blood pressure, which suddenly dropped in the standing position, hence, use 
of infusion of ephedra very was considered effective in this condition in 
which the person's blood pressure became normal but it can not be used 
for a long time because it may also increase blood pressure [2, 34]. 


Ephedra has catechols constituents, which are alpha-receptor activity 
(vaso-constriction), due to which it can have vasoconstriction effects. 


3.2. Respiratory System 


Ephedra has two basic constituents’ ephedrine and pseudoephedrine. 
Both of them gives a very good effect on the respiratory system. The 
constituents help to widen the lung's airway and also help to cure the nasal 
decongestant. 


Additionally, they also help to treat asthma because these constituents 
have b-adrenergic receptors, which relax airway smooth muscles, hence, 
ephedra is noticed to be effective in asthma cases [2]. 


Ephedra has been effective in the respiratory system, especially in 
asthma that is why it is also known as DAMA BOTI. 


3.3. Anti-Inflammatory 


Ephedra also helped to reduce swelling. It has ephedranine A & B 
chemicals, which contain anti-inflammatory effects and reduce swelling, 
especially when applied in joint pain conditions like Rheumatoid Arthritis 
and Osteoarthritis. Paste of bark of ephedra should be applied on the 
affected area, which may give relief in pain and also in swelling [3]. 


Due to its effect on the joint, it is called Joint Pine. 


A 
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3.4. Anti-Tumor Activity 


Ephedra also has anti-tumor activity. It has two chemicals, namely, 
lignans, flavones both of which have anti-oxidant and anti-tumor properties, 
which suppress the growth of the tumor and also act as anti-bacterial 
medicine [2]. 


3.5. Renal System 
Ephedra has another special property, which is Diuretic [1, 2]. 


It chemical ephedrine has multiple actions and one of them is Diuretic 
action. It produces a positive effect on renal system diseases such as renal 
calculi, urethritis, and UTI. , and those conditions in which dieresis is 
necessary through the kidney [1, 34] 


3.6. Side Effects 


Normally, any drug either it is herbal or synthetic, its excessive use can 
be harmful. 


Prolong use of joint pine could elevate blood pressure because it has 
vasoconstriction property, which also helps to reduce lower blood pressure 
[2, 3]. 


It contraindicated in high blood pressure. 


Some noticeable side effects include, vomiting, nausea, and minor 
abdominal pain [1, 35]. 


Some researches is also used in Obesity but in my point of view it not 
be useful for slimming medicine because it might be caused hypertension 
[1, 35]. 


4. CONCLUSION 


Traditional medicine has employed E. gerardiana for its 
bronchodilatory, diaphoretic, diuretic, and thermogenic properties. The 
bronchodilatory effects of ephedrine have been extensively studied, leading 
to its inclusion in various over-the-counter medications for respiratory 
conditions. However, it is important to note that ephedrine-containing 
products have been associated with potential side effects and safety 
concerns. Adverse effects such as increased heart rate, elevated blood 
pressure, and central nervous system stimulation have been reported. 
Therefore, caution should be exercised when using E. gerardiana or 
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products containing ephedrine and it is advisable to consult with healthcare 
professionals before its usage. Further research is needed to elucidate the 
precise mechanisms of action of Ephedra gerardiana's bioactive 
compounds, optimize dosage regimens, and explore potential drug 
interactions. These studies would contribute to ensuring the safe and 
effective use of E. gerardiana in clinical settings and enhancing our 
understanding of its therapeutic potential. Concluding, E. gerardiana is a 
valuable medicinal plant but it should be used carefully, considering the 
potential risks associated with ephedrine-containing products. 
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